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ABSTRACT 

When combined with infrared observations with the Spitzer telescope (3 to 160 yum), the 
Herschel Space Observatory now fully samples the thermal dust emission up to 500 yum and 
enables us to better estimate the total infrared- submm energy budget (Ltir) of nearby galax- 
ies. We present new empirical calibrations to estimate resolved and integrated total infrared 
luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. 
We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed 



with Herschel. We perform a resolved SED modelling of these objects using the Draine & Li 



( |2007| dust models and investigate the influence of the addition of SPIRE measurements in 



the estimation of Ltir. We find that using data up to 250 yum leads to local Ltir values con- 
sistent with those obtained with a complete coverage (up to 500 jim) within +10% for most 
of our resolved elements. We then study the distribution of energy in the resolved SEDs of 
our galaxies. The bulk of energy (30-50%) is contained in the [70-160 yum] band. The [24-70 
Jim] fraction decreases with increasing metallicity. The [160-1100 yum] submillimeter band 
can account for up to 25% of the Ltir in metal-rich galaxies. We investigate the correlation 
between TIR surface brightnesses/luminosities and monochromatic Spitzer and Herschel sur- 
face brightnesses/luminosities. The three PACS bands can be used as reliable monochromatic 
estimators of the Ltir, the 100 yum band being the most reliable monochromatic tracer. There 
is also a strong correlation between the SPIRE 250 yum and Ltir, although with more scatter 
than for the PACS relations. We also study the ability of our monochromatic relations to repro- 
duce integrated Ltir of nearby galaxies as well as Ltir of z~l-3 sources. Finally, we provide 
calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from 
a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated 
uncertainties. 
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1 INTRODUCTION 

Interstellar dust obscures our view of the star formation sites in 
galaxies. Indeed, 30 to 50% of the starlight emission is thermally 
reprocessed by dust, and re-emitted at infrared (IR) wavelength 
( Draine 2003| |Tielens|[2005| ). This wavelength regime enables us 
to directly investigate the dust physics and indirectly probe the star 
formation activity obscured by dust within galaxies and is thus cru- 
cial to understand how galaxies evolve and recycle their interstel- 
lar material. The total bolometric IR emission Ljir constitutes the 
emission of all the dust-enshrouded stellar populations (but can also 
include emission from Active Galaxy Nucleus or AGN) and is one 
of most reliable tracers of the star formation obscured by dust. Sev- 
eral studies have thus derived calibrations of the star formation rate 
(SFR) based on Ltir ( Kennicutt 1998, Perez-Gon zalez et al.|2006| 
IKennicutt et al.|2009, Kennicutt & Evans 2012J . 

Ltir can be estimated by combining multi- wavelength obser- 
vations sampling the thermal dust emission from mid-IR to submil- 
limeter wavelengths and integrating the emission directly or using 
realistic dust models to interpolate the data. Unfortunately, many 
galaxies do not benefit from a complete sampling of their Spec- 
tral Energy Distribution (SEDs), which prevents the modelling of 
their SEDs and thus a correct estimate of their Ljir. Previous works 
have thus provided calibrations of the Ltir using monochromatic 
IR wave bands or a combination of IR wave bands. For instance, 
[Sanders & Mirabel] (T996l l or |Sanders et al.| ( |200 3) provided a re- 
lation to derive the Ltir of luminous IR galaxies using the 4 IRAS 
(Infrared Astronomical SateUite) filters at 12, 25, 60 and 100 jim. 
[Dale & Helou| ( [2002 | updated this relation using a combination of 
the three Spitzer/MIFS wavelengths (24, 70 and 160 //m), match- 
ing their modelled Ltir with very good accuracy. More recently, 
[Boquien et al. (2010j used the |Dale & Helou ( 2002 ) relation to es- 
timate the Ltir from other Spitzer bands, the 8 and 24 fim bands in 
particular. 

The good resolution of the two IR-submillimeter instruments 
PACS and SPIRE onboard the Herschel Space Observatory opens 
a new window on how to quantify Ltir at local scale. The reso- 
lutions of the SPIRE instrument match that of the MIPS instru- 
ment of the Spitzer Space telescope, the predecessor of Herschel. 
Indeed, the resolution of SPIRE at 250 fim (~18'') is similar to that 
of SpitzerlMWS 70 jim and the resolution at 500 jim (~36'') is sim- 
ilar to that of MIPS 160 fim. Furthermore, Herschel data enable a 
more complete coverage of the peak of the thermal dust emission 
and of the submm slope of nearby galaxies up to 500 jim, allowing 
a refinement of our measurements of the Ltir. 

Combining Spitzer and Herschel bands, [Boquien et aL] ( |2011| ) 
derived resolved estimators of the total infrared brightness in the 
galaxy M33. In this study, we aim to similarly model resolved Ltir 
for a wider sample of galaxies, investigate the relation between Ltir 
and diff'erent Herschel bands, study how those relations evolve with 
the galaxy characteristics as well as provide recipes to obtain reli- 
able Ltir predictions from a large choice of single or combined 
wavelengths. As previously mentioned, many authors have stud- 
ied relations between IR monochromatic fluxes and Ltir as well. 
Recipes like those provided by [Sanders & Mirabell ( [1996| l and jDalej 
[& Helouj ( |2002j ) are often used in the literature as a proxy for the 
derivation of Ltir. This estimated Ltir is then used to derive cal- 
ibrations from monochromatic luminosities ( [Boquien et al.|2010| 
[Elbaz et al.|2011[ among others). One of the advantages of the ap- 
proach we follow in this study is that we now have access to the 
whole coverage of the thermal dust IR emission in our nearby ob- 
jects with Herschel. Our Ltir will be directly modelled using the 



IR observations and a realistic dust SED model, which limits un- 
certainties and biases linked with previous calibrations. 

We perform this analysis using the Herschel data of ~60 
nearby galaxies observed as part of the KINGFISH (Key Insights 
on Nearby Galaxies: A Far-Infrared Survey with Herschel; |Kenni-[ 
[cutt et al.[2011j ) programme. The paper is organised as follows. We 
present the sample, Spitzer and Herschel data in Section 2. As we 
aim to derive resolved estimators of Ltir, we would like to use the 
highest resolution available. The first step of this study is thus to 
determine the best compromise between resolution and sufficient 
constraint on the Ltir estimates (Section 3). We then analyse the 
distribution of the total infrared energy with wavelength on a lo- 
cal basis in Section 4. We present a calibration of the TIR surface 
brightnesses/luminosities using single Spitzer or Herschel bands 
in Section 5 as well as calibrations from a combination of vari- 
ous bands in Section 6. Because the metallicity and the hardness 
of the radiation field are parameters that strongly aff'ect the far-IR 
emission and the range in dust temperature from galaxy to galaxy, 
throughout the paper, we investigate how our relations and their 
reliability evolve with global or local galaxy properties. 



2 A MULTI- WAVELENGTH MAPPING 
2.1 The sample 

We obtain the Herschel data (PACS and SPIRE maps) as part of 
the Herschel key programme KINGFISH. This sample provides a 
unique opportunity to study the relation between Herschel bands 
and total IR luminosities. The sample comprises 61 galaxies, prob- 
ing a wide range of galaxy types (from elliptical to irregular galax- 
ies) and various star-formation activities, from active star-forming 
regions to more q uiescent ISM, with global SFRs ranging from 
10-3 to 7 Mo yr-i ( jHowell et al1 ( [201 0') even estimate a SFR of 23 



Mq yv~^ for the luminous IR galaxy NGC 2146). It also includes 
galaxies hosting low-luminosity AGN. The KINGFISH galaxies 
are located between 3 and 31 Mpc, leading to ISM resolution ele- 
ments of 0.2 to 2.6 kpc at the resolution of SPIRE 250 jim (FWHM 
of the PSF: 18''), the resolution at which we work in the follow- 
ing study (see Section 3 for justification). The KINGFISH galaxies 
also probe various metallicities. We use the metallicities tabulated 
by [Kennicutt et al.| ( |201 Ij ) who provide two metallicities per galaxy, 
one derived from the theoretical calibration of [Kobu lnicky & Kew- 
ley ( 2004 ), the other from the empirical calibration of Pilyugin & 
,Thuan (2005 K Here we use the latter; with this calibration, oxygen 
abundances (defined as 12+log(0/H)) range from 7.54 for the low- 
metalHcity galaxy DDO 154 to 8.9 for the galaxy NGC 3077. We 
note that metallicity gradients are observed in some of the KING- 
FISH galaxies (Moustakas et al. 2010 ) but only few gradients are 
currently well constrained. In this paper, we adopt the same metal- 
licity in each resolution element of a given galaxy, equal to that de- 
termined globally for the galaxy using the |Pilyugin & Thuan| ( [2005| ) 
calibration. 

Previous and on-going studies are also analysing global 
and local SED models of galaxies of the KINGFISH sample 
using the Herschel data. [Dale et aT] p012| ) present the Herschel 
far-IR and sub-millimeter photometry of the KINGFISH survey 
as well as integrated SED models of these objects from which 
total dust masses are in particular derived. [Skibba et al.[ p011| ) 
also compare the global emission from dust and from stars in the 
same sample and analyse how the dust- to- stellar flux ratio varies 
with properties such as morphology, Ltir or metallicity. Using 
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local modified blackbody models for a sample of the KINGFISH 
galaxies, [Galametz et al.| ( [20T2| investigate the physical properties 
(temperature, emissivity) of the cold dust phase and associated 
uncertainties. Finally, Aniano et al.| ( [201 2 ) present a pixel-by-pixel 
SED modelling and a mapping of the dust and radiation field 
properties for the two spirals NGC 628 and NGC 6946 using the 
[Draine & Li| ( |2007| ) dust models. This local modelling will be 
extended to the whole KINGFISH sample in Aniano et al. (in prep). 



2.2 Herschel maps 

The Herschel/FACS instrument ( |Poglitsch et al.|[2"010| ) provides 
maps with FWHMs of the PSFs of 5776x5? 46, 6? 69x6789 and 
12713xl(y.'65 at 70, 100 and 160 fim respectively for the chosen 
scan speed (20'7s). Observations of the KINGFISH galaxies with 
this instrument were obtained with 15' long cross-scans (perpen- 
dicular scans). From raw data to Level 1, the processing of PACS 
data follows the main steps of the recommended standard proce- 
dure for steps of pointing association, conversion to physical units 
or flat-fielding. Glitches are removed using a second-level deglitch- 
ing method based on a comparison of individual readouts with a 
reference sky value at the same position that allows us to detect 
outlier values. We refer to [Kennicutt et al. | ( |20 ll\ for more details 
on the initial processing of the data within the Herschel Interactive 
Processing Environment (HIPE, version 8). 

We use the Scanamorphos technique (version 16.9) to process 
the data from these Level 1 data, correct for 1/f noise and project 
the pixel timelines in the sky in order to build the final maps. 
Scanamorphos in particular subtracts the brightness drifts caused 
by the low-frequency noise using the redundancy built in the obser- 
vations ( |Roussel|2012| ). The final pixel sizes of our PACS data are 
174, 177 and 2785 at 70, 100 and 160 jim respectively. The PACS 
calibration uncertainties are ~5 9^ 

The Herschel/SFIRE instrument ( [Griffin et ari[20T0l ) pro- 
vides maps with FWHMs of the PSFs 1873x17'', 2477x2372 and 
37"x3374 at 250, 350 and 500 yum respectively. Observations were 
obtained in scan mode. The data reduction was performed from 
raw data with the HIPE environment. We refer to the KINGFISH 
Data Products Delivery User's Guid^for details on the data reduc- 
tion. The SPIRE maps used in this study are built with a nearest- 
neighbor projection on sky and averaging of the time ordered data. 
The final pixel sizes of our SPIRE data are 6", 10" and 14" at 
250, 350 and 500 yt/m respectively. Calibration uncertainties are es- 
timated to be ~7% for the three wave band^ 

We refer to|Kennicutt et al.| ( [20n] ), |Engelbracht et aLlpOTO] ) 
and [Sandstrom et al.| ( |20I0| )"for more details on the KINGFISH 
sample, the observation strategy and the different steps of the data 
processing. We do not include the KINGFISH galaxies DDO 154, 
DDO 165, Holmberg I and NGC 1404, since they are barely de- 
tected with Herschel (upper limits in the global flux catalogue of 
[Dale et al.||2012| ). We also note that PACS observations for the 
galaxy NGC 584 are contaminated by emission from Jupiter. This 
galaxy is thus also excluded from the following analysis. 



^ http://herschel.esac.esa.int/Docs/PACS/html/pacs_om.html 

^ http://herschel.esac.esa.int/UserReducedData.shtml 

^ http://herschel.esac.esa.int/Docs/SPIRE/html/spire_om.html 



2.3 Spitzer maps 

Most of the KINGFISH galaxies have been observed with 
Spitzer/IRAC and MIPS as part of the SINGS programme (Spitzer 
Infrared Nearby Galaxies Survey; ' Kennicutt et al.|[2003| ). Four 
galaxies of the sample are drawn from other Spitzer surveys: 
IC 342, NGC 5457 (MlOl), NGC 2146 and NGC 3077. IRAC ob- 
serves at 3.6, 4.5, 5.8 and 8 jim with PSF FWHMs of V!7, V!7, 
V!9 and 2" respectively. The IRAC images are reduced using the 
SINGS Fifth Data Delivery pipelin^ Maps are multiphed by 0.91, 
0.94, 0.66 and 0.74 at 3.6, 4.5, 5.8 and 8 yt/m respectively to account 
for extended- source flux calibratiorj^ MIPS observes at 24, 70 and 
160 jim with FWHMs of the PSFs of 6'^ 18'' and 40'' respectively. 
Because of their lower resolution compared to PACS 160 yt/m maps, 
we do not use the MIPS 160 //m maps in the following study. Galax- 
ies of the LVL (Local Volume Legacy) survey are reduced using the 
LVL pipelin^ Galaxies that are not part of the LVL survey were 
re-processed using the LVL reduction technique for consistency. 



3 INFLUENCE OF SPIRE DATA ON THE Ljir MAPS 

We aim to derive Lrm estimators using Spitzer and Herschel bands 
(as monochromatic or combined tracers). To make the most of 
the good resolution of Herschel and work at the highest resolu- 
tion available, we first investigate in this section how SPIRE wave- 
lengths influence the determination of the hjiR. We thus derive Ljir 
maps using data between 3.6 and 160 jum, between 3.6 and 250 /um 
or between 3.6 and 500 jum to quantify the diff'erence in the global 
and local Ljir values. 

3.1 Obtaining the I^tir maps 

In a first step, we subtract the sky background of Spitzer and Her- 
schel maps using a plane- subtraction technique (see [Aniano et al.| 
|2012| and Aniano et al., in prep). We then use the convolution ker- 
nels provided by [Aniano et al.| ( |2011| ) to convolve the Spitzer and 
Herschel maps to three diff'erent resolutions: 

- PACS 160 jjim resolution : We convolve the IRAC 3.6 /im, 
4.5 yL/m, 5.8 yL/m and 8.0 yt/m, the MIPS 24 yt/m, the PACS 70 yt/m and 
PACS 100 yL/m maps to the resolution of PACS 160 yt/m (FWHM ~ 
12") and regrid them to a common pixel size of 4" (original pixel 
size of the PACS 160 fim image), 

- SPIRE 250 jjim resolution : We convolve the IRAC 3.6 yum, 
4.5 yL/m, 5.8 yL/m and 8.0 yL/m, the MIPS 24 yL/m, the PACS 70 yL/m, 100 
yL^m and 160 yL/m maps to the resolution of SPIRE 250 jjm (FWHM 
~ 18") and regrid them to a common pixel size of 6" (original pixel 
size of the SPIRE 250 yL/m image), 

- SPIRE 500 Jim resolution : We convolve the IRAC 3.6 yL/m, 
4.5 yL/m, 5.8 yL/m and 8.0 yL/m, the MIPS 24 and 70 yL/m, the PACS 70 
yL/m, 100 yL/m and 160 yL/m, the SPIRE 250 and 350 yL/m maps to the 
resolution of SPIRE 500 jim (FWHM ~ 36") and regrid them to 
a common pixel size of 14" (original pixel size of the SPIRE 500 
yL/m image). 

Using our maps convolved at three diff'erent resolutions - and 

^ http://data.spitzer.caltech.edu/popular/sings/20070410_enhanced_vl/Documents/ 
sings_fifth_delivery_v2.pdf 

^ http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/iracinstrumenthandbook/ 
6 http://irsa.ipac.caltech.edu/data/SPITZER/LVL/LVL_DR5_v5.pdf 
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consequently providing three different coverages - we perform lo- 
cal SED fits using the |Draine &Ti| < [2007t (hereafter [DL07]) dust 
models in order to match the observed fluxes in each resolution 
element. We refer to |Aniano et al.| { [2012| ) for a full description 
of the pre-data treatment (convolution or background subtraction 
steps, production of the uncertainty maps for each bands etc.) and 
the resolved SED modelling process (description of the model, as- 
sumption on parameters etc.). The SED modeling uses Spitzer and 
Herschel bandpasses directly and alleviates the need for colour- 
corrections. In most regions, stellar emission dominates the emis- 
sion of the two first IRAC bands (3.6 and 4.5 jim) while dust dom- 
inates observations beyond 4.5 jim and at least up to 500 jim (see 
[Engelbracht et al.|20 08 for instance). To account for thermal dust 
emission only, we subtract the contribution of stellar emission to 
the short wavelengths during the modelling process. We approx- 
imate the stellar emission at /I > 3 yt/m by scaling a blackbody 
function, using a representative photospheric temperature of 5000K 
( |Bendo et al.|2006[[Drame et al.|2007l ). 

l^TiR measures the total dust emission and is obtained by in- 
tegrating the SED in a v-Ly space. In this paper, we define I^tir 
as: 



Table 1. Integrated I^tir values for different SED coverages. 



^1100 
*J3 fim 



1 100 lum 



Ly dv 



(1) 



We thus integrate the SEDs from 3 to 1100 jim to obtain 
the Ltir in each resolved element. We note that our resolved 
SED models have a logarithmically- spaced wavelength grid 
of ~350 values from 3 to 1100 yt/m. We use the IDL function 
INT_TABULATED (5-point Newton-Cotes formula) to perform the 
integration. We restrict ourselves to resolved elements that 1 ) are 
located within the elliptical apertures used by [Dale et al.| ( |2012] ) to 
perform the global photometry, 2j do not contain contamination 
from foreground stars or known background galaxies along the 
line-of-sight and 3) have a 3-cr detection in all the bands used for 
the modelling. This leads to Ljir maps of our KINGFISH galaxies 
obtained at three different resolutions (and for three different SED 
coverages). 

In the rest of this paper, we use the nomenclature: 

- Ltir pi 60 • the Ltir modelled with data constraining the SED 
from 3.6 to 160 fim (at PACS 160 yum resolution), 

- Ltir 5250 ■ the Ltir modelled with data constraining the SED 
from 3.6 to 250 fim (at SPIRE 250 yt/m resolution), 

- Ltir 5500 • the Ltir modelled with data constraining the SED 
from 3.6 to 500 yt/m (at SPIRE 500 yt/m resolution). We consider 
our Ltir 5500 maps as our "reference" maps, since resolved SEDs 
were modelled with the most complete coverage of the thermal 
dust emission. 

Submm excess or AGN contribution - Ground-based data and 
Herschel observations at submm wavelengths have helped us to 
better investigate the properties of the coldest phases of dust. A 
flattening of the submm slope or an excess compared to fits per- 
formed without submm data is very often detected in metal-poor 
galaxies ( GaUi ano et al. 2003||20051|Marleau et al.l2006[[Gaiametz| 
|etar.|2009t|O'Halloran et al.f 2010, Bot et al.'201ol among others). 
Various hypotheses have been investigated to explain this excess: 
emission from a shielded cold dust reservoir UGalHano et al.|2005[ 
[Gal ametz et al.'2009), temperature-emissivity dependence of dust 
grains (Meny et al. 2007 , Paradis et al. 2010 ), "spinning dust" emis- 
sion |Murphy et al. [20 10[ [Planck Collaboration|2011j| ) and recently 
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DDO 053 


1.24x10^ 


8.83 -^-^xlO^ 


1.08 


1.04 


Holmberg II 


7.08x10^^ 


6.09 -^•^xlO'^ 


0.98 


0.96 


IC 342 




1.47 -^-^xlO^^ 


1.06 


1.03 


IC 2574 


2.04x10^ 


1.69 -^-^xlO^ 


0.98 


0.99 


M81 dwB 




3.54 -^-^xlO^ 


1.02 


1.00 


NGC 0337 


1.25x10^^ 


1.06 -^-^xlO^^ 


1.04 


1.06 


NGC 0628 


7.78x10^ 


6.89 -^-^xlO^ 


1.04 


1.02 


NGC 0855 


3.91x10^ 


3.55 ±0-2x10^ 


1.02 


1.01 


NGC 0925 


4.18x10^ 


3.68 ^O'^xlO^ 


0.97 


0.98 


NGC 1097 


4.56x10^^ 


4.16-O-^xlO^o 


0.99 


1.01 


NGC 1266 


2.63x10^^ 


2.45 ^^•'^xlO^^ 


1.04 


0.98 


NGC 1291 


2.76x10^ 


2.70 -^-^xlO^ 


0.93 


0.97 


NGC 1316 


6.83x10^ 


5.95 ±0-3xio9 


1.08 


1.02 


NGC 1377 




1.26 ±01x10^0 


1.02 


1.00 


NGC 1482 


5.05x10^^ 


4.53 ±0-9x1010 


0.97 


1.03 


NGC 1512 


3.51x10^ 


3.43 -^-^xlO^ 


1.05 


1.04 


NGC 2146 




1.30 ^O-^xlO^i 


0.98 


0.95 


NGC 2798 


4.15x10^^ 


3.36 ^O-^xlO^o 


0.97 


1.02 


NGC 2841 


1.11x10^^ 


9.22 -^-^xlO^ 


1.04 


1.03 


NGC 2915 


4.19x10^^ 


3.36 ^O'^xlO^ 


1.14 


1.05 


NGC 2976 


8.46x10^ 


7.57 ±0-7xio8 


1.00 


1.00 


NGC 3049 


3.65x10^ 


3.21 ^O'^xlO^ 


1.06 


1.08 


NGC 3077 




7.31 ±0-^x10^ 


0.98 


1.02 


NGC 3184 


1.07x10^^ 


8.45 -^-^xlO^ 


1.03 


1.02 


NGC 3190 


6.61x10^ 


5.96±iOxio9 


1.05 


1.02 


NGC 3198 


8.83x10^ 


7.27 -^-^xlO^ 


1.02 


1.02 


NGC 3265 


2.83x10^ 


2.49 -^-^xlO^ 


1.03 


1.02 


NGC 3351 


7.81x10^ 


6.86±iOxio9 


1.03 


1.03 


NGC 3521 


3.24x10^^ 


3.15 ±0-2x1010 


1.04 


1.03 


NGC 3621 


7.90x10^ 


6.80 ±0-^x10^ 


1.03 


1.01 


NGC 3627 


2.66x10^^ 


2.53 ±0-3x1010 


1.00 


1.01 


NGC 3773 


6.69x10^ 


5.04 ±1-^x10^ 


1.07 


1.06 


NGC 3938 


1.87x10^^ 


1.57 ±0-2x1010 


1.02 


1.02 


NGC 4236 


4.68x10^ 


3.85 ±0-^x10^ 


1.05 


1.04 


NGC 4254 


4.21x10^^ 


3.57 ±0-2x1010 


1.03 


1.02 


NGC 4321 


3.38x10^^ 


2.88 ±0-2x1010 


1.02 


1.01 


NGC 4536 


2.12x10^^ 


2.06 ±0-2x1010 


0.98 


1.05 


NGC 4559 


3.01x10^ 


2.51 ±0-3x10'^ 


0.97 


1.00 


NGC 4569 


1.41x10^^ 


1.37 ±0-3x1010 


1.05 


1.04 


NGC 4579 


1.20x10^^ 


9.22 ±0-^x10^ 


1.06 


1.06 


NGC 4594 


3.42x10^ 


3.03 ±0-3x10^ 


1.03 


1.03 


NGC 4625 


6.19x10^ 


5.20 ±0-^x10^ 


1.00 


1.01 


NGC 4631 


2.54x10^^ 


2.04 ±0-2x1010 


0.98 
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Figure 1. Integrated I^tir sings from |Draine et al.|{2007} (top), Ltir pm (middle) or Ltir 5250 (bottom) compared to our reference Ltir ssoo- Integrated 
Ltir are expressed in Lq. For the top panel, we indicate with filled circles the luminosities Ltir sings derived including SCUBA data at 850 fim in the SED 
fitting and with empty circles the galaxies for which SCUBA data were not available. 



magnetic dipole radiation from magnetic nanoparticles ( [Draine &] 
|Hensley||2012| ). In the following study, we do not include data at 
wavelengths greater than 500 fim (wavelength at which submm ex- 
cess starts to be detected) but we include modelled emission out to 
11 00 yL/m in our estimations. Our estimations of Ltir would not take 
the submm excess into account, if any. Nevertheless, as quantified 
further in this paper (Section 4), we expect thermal dust emission 
above 500 fim to be negligible in the bolometric infrared energy 
budget of our galaxies, especially for low-metallicity environments 
where the excess is usually detected. In addition, a large number 
of the KINGFISH galaxies show nuclear emission indicating ex- 
citation by a non stellar continuum but no dominant AGN (except 
NGC 1316). We do not expect the AGN contribution to signifi- 
cantly afifect our estimates of Ltir. 



3.2 Comparison of the integrated Ltir 

Ltir being a linear quantity, we do not expect strong dififerences be- 
tween integrated Ltir obtained from a SED model fitting the galaxy 
as one big pixel or the Ltir we obtained by summing the resolved 
Ltir. To check this assumption, we derive the integrated Ltir using 
the two techniques. We find that for each galaxy, the two integrated 
Ltir values diflfer by 5% at most whatever the wavelength cover- 
age, which is comparable to the error bars on these quantities. For 
the rest of the paper, integrated Ltir values are thus obtained by 
summing the resolved Ltir. 

We want to investigate how integrated Ltir values vary when 
we include SPIRE data in the SED modelling. We thus derive inte- 
grated Ltir by summing the resolved Ltir in the |Dale et al.H2012] l 
apertures obtained using a [3-160yL/m] coverage, a [3-250yLim] cov- 
erage and a [3-500//m]. Table [T] lists the integrated Ltir values 
of our sample for the three diflferent coverages. Global errors are 
obtained by summing the resolved uncertainties estimated during 



Galametz et al, 
NGC 628 



NGC 6946 



+ i5°5r 

+ 15°48' 
+ 15°45' 
+ 15°42' 



+ 15°48' 
+ 15°45' 
+ 15°42' 

+ 15°5r 
+ 15°48' 
+ 15°45' 
+ 15°42' 




LtIR P160 conv S500 ^ 



if * 



(Lpieo Lggoo) 



0.3 



I 



^S500^ 



-0.3 



LtIR S250 conv S500 ^ 




-(Ls250 Lg5oo) 



0.3 



-0.3 : 



^S500 



f60°12' 



f60°6' 



f60°3' 



+ 60°12' 



f60°9' 



f60°6' 



f60°9' 



f60°6' 



f60°3' 



^1^1 H^i^ i 



7 Lt 



7 



LtIR P160 conv S500 




I 



I 



7 L, 



'TIR S250 conv S500 




^ - L^iR S500 

- 20''35"'12' 20''34"'24' 



20''35"'12' 20"34'^24' 20^35"'12' 20''34"'24' 



Figure 2. Comparison of I^tir maps of NGC 628 and NGC 6946 for different resolutions. For each galaxy: Top: I^tir p\60^ ^tir 5250 and hjiR ssoo maps 
(Lq kpc^, log scale). Middle: hjjR pi^q and hjiR 5250 maps convolved to the SPIRE 500 /urn resolution (L© kpc^, log scale). The last column is unchanged. 
Bottom: Relative difference between the convolved Ltir pi60 or Ltir 5250 maps and the reference Ltir 5500 map. 
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Figure 3. Resolved differences between 1^tirpi6Q or ^tirsisq and the reference I^tirssqq as a function of the TIR brightness Stirssqq- Galaxies are 
metallicity coded: purple for metal-rich, red for low-metallicity (colour scale in terms of 12-i-log(0/H) in the top panel). For each galaxy, we sort and average 
the data 3-by-3 for clarity. The black points indicate the median of the complete distribution per TIR brightness bins with 1-cr error bars. A histogram of the 
relative difference for the full sample of pixels is shown on the right of each plot. 



the SED modelling process, with a median offset of ~12% for 
the integrated I^tir ssqq- We normalise the integrated I^tir pieo and 
Lr/i? 5250 to l^TiR S500 for comparison. As part of the SINGS project, 
[Draine et ah] (2007 ) modelled most of the KINGFISH galaxies us- 
ing Spitzer fluxes (thus an SED coverage up to 160 jim) and the 
[DL07] dust models and derive global dust luminosities (hereafter 
^TiR sings)- They also include SCUBA fluxes at 850 fim, when 
available, for a small subsample of their objects. In this paper, 
we use the distances provided by Kennicutt et al.| ( |2011), some of 
them being different from those used in the study of [Praine et al.| 
([2007 ). We thus first rescale the hnR sings values to the distances 
we choose to use and add these corrected values to Table[T]for com- 
parison. 

Figure [T] illustrates how the integrated Lrm sings (top), 
^TiR PI 60 (middle) or Ltir siso (bottom) compare with the inte- 
grated hjiR obtained in our complete coverage case Ltir ssoo- For 
the top panel, filled circles indicate when SCUBA data at 850 
fim were used in the fit to determine hnR sings ^ empty circles 
when SCUBA data were not available. Our integrated Ljir ssoo 
estimates are close to the Ltir sings estimates within 10-15%, even 
if systematically lower than the values derived from Spitzer data 



only. Comparing the top and middle panel, we observe differences 
(a shift) between integrated Ltir sings and Ltir p\60^ both obtained 
using the same coverage up to 160 yum. As previously mentioned, 
resolution effects are not sufficient to explain such a shift between 
SINGS values and our values. The difference is thus probably 
related to the use of PACS data in the fit compared to the MIPS 
data used previously (see [Aniano et al.||2012l for discussion on 
PACS-MIPS photometry disagreement). Integrated Ltir values 
obtained with our three different coverages (Ltir pieo^ Ltir 5250. 
Ltir 5500) are very similar. The agreement is within 8% between 
the integrated Ltir pm or Ltir siso and the reference Ltir 5500 
(except for NGC 2915 for which Ltir pieo is higher by 14%). 



3.3 Comparison of the resolved Ltir 

We now probe the variations in the Ltir estimates driven by the 
different coverages on a resolved scale. We convolve the Ltir pm 
and Ltir 5250 maps of the KINGFISH galaxies to the resolution of 
SPIRE 500 yL/m and compare them to our reference Ltir ssoo map. 
Figure |2] gives an example of this comparison for the two spiral 
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galaxies NGC 628 and NGC 6946. The top panels show the Ltir 
maps at the original resolution (PACS 160, SPIRE 250 and SPIRE 
500 yL/m from left to right) in L© kpc~^. The middle panels show 
these maps convolved to the SPIRE 500 fim resolution (the last 
column is unchanged). The bottom panels finally compare the con- 
volved hjiR pieo or p250 maps with the reference Ljir ssoo map. 
We remind the reader that we restrict our study to resolved ele- 
ments with a 3-cr detection in PACS 160 and SPIRE 250 and 500 
fim bands. We can still distinguish the structure of NGC 628 and 
NGC 6946 on the bottom left panels that show the diff'erence be- 
tween the hjiR P160 maps and the reference maps Ljir ssoo- In both 
galaxies, we observe that the absence of submm (SPIRE) data leads 
to an underestimation of the Ltir in bright regions (blue structures 
in Figure |2]) of up to 10-15%, and an overestimation in the outer 
part of the galaxies, so at low-surface brightnesses (red structures), 
of up to 30%. The Ltir maps obtained with data up to SPIRE 250 
and with data up to SPIRE 500 yt/m are, on the contrary, very sim- 
ilar, as illustrated by the small diff'erence residuals of those maps 
(bottom middle panels). 

Figure [3] gathers, for the complete sample, the resolved rela- 
tive diff'erences between Ltir pm (top) or Ltir 5250 (bottom) with 
the reference Ltir 5500 as a function of the TIR surface brightness 
^TiR 5500 (in W kpc~^, log scale). Galaxies are metallicity-coded, 
from low-metallicity in red to high-metallicity in dark purple 
colour. Discrepancies between resolved Ltir pm and Ltir 5500 can 
be significant as illustrated in Fig. [3] (top). Ltir is under-predicted 
when modelled with data up to PACS 160 fim in high surface 
brightness regions (drop of the median for log Stir 5500 > 34.5 
W kpc~^), and over-predicted in low-surface brightness regions, 
similar to the trends observed for NGC 628 and NGC 6946 
(Fig. |2]). Using data up to 250 yt/m resolution, we can observe 
a slight under-prediction of Ltir at 33 < log Stir 5500 < 34 (in 
W kpc~^) and an over-prediction at log Stir ssoo > 34 but Ltir 5250 
and Ltir ssoo are nevertheless close (within +10% for most of our 
resolved elements). For both resolution, determining the Ltir with 
accuracy seems difficult at very low Stir (< 33 W kpc~^), due to 
the large uncertainties on the flux measurements in those regions. 
In conclusion, using data up to 250 fim seems to be the best 
compromise between sufficiently constraining the submm slope in 
order to obtain a Ltir consistent with that obtained with a complete 
coverage of the dust thermal emission, while still keeping a good 
working resolution. 

We choose to work at the resolution of SPIRE 250 /dm for the 
rest of this study. At this resolution, resolved Ltir have uncertain- 
ties of 10-15% on average in the resolved elements we select. In the 
Appendix (Figure [dT] ), we show the surface brightness Stir maps 
obtained at SPIRE 250 fim resolution for the full KINGFISH sam- 
ple. Maps are in L© kpc~^ (log scale). We remind the reader that 
the FWHM of the SPIRE 250 PSF is ~18'^ Our final maps have a 
pixel size of 6", which corresponds to ISM elements of 88 pc for 
the closest galaxy of the sample Holmberg II (3.05 Mpc) and 890 
pc for the furthest galaxy of the sample NGC 1266 (30.6 Mpc). 



4 THE INFRARED TO SUBMM DISTRIBUTION OF 
LUMINOSITIES 

For normal star-forming galaxies, |Dale & Helou| p002| ) studied 
the distribution of the infrared energy budget, namely how much 
energy emerges in various wavebands, and investigate its depen- 
dence on the star formation activity. [Dale et al.|p009| ) also present 
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Figure 4. Resolved fraction of Ltir emitted in different wavebands ([3-24], 
[24-70], [70-160] and [160-1100]) as a function of the fv(70 /im)/fv(100 
fim) far-IR colour. Galaxies are metallicity coded: purple for metal-rich, red 
for metallicity-poor (colour scale in the top panel in terms of 12-i-log(0/H)). 
For each galaxy, we sort and average the resolved element 5-by-5 for clarity. 
Global luminosity fractions are overlaid with black circles. A log -linear fit 
to the data in the bottom panel is shown by the grey dashed line. 



monochromatic-to-bolometric ratios and analyse their dependence 
with fv(70 jum) / fy (160 /um) ratios and morphologies. In a simi- 
lar fashion, we analyse the distribution of the total infrared energy 
with wavelength for the KINGFISH sample, now using our Her- 
schel data on a resolved scale. We derive the fractions of the Ltir 
emitted in 4 different wavebands ([3-24 yum], [24-70 yum], [70-160 
yL/m] and [160-1100 /urn]): we integrate the local SED models over 
the stated limits and divide these values by the resolved Ltir. 

Figure [4] displays these fractions as a function of the 70-to- 
100 fim flux density ratio (hereafter 70/100 colour), which pa- 
rameterises various star formation activities and is also a good 
proxy for dust temperatures. This ratio also correlates well with 
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the starlight intensity. Galaxies are metallicity coded, from dark 
purple for metal-rich to red for metal-poor. We overlay integrated 
values for comparison. We obtain averaged 70/100 colours (x-axis) 
by summing the 70 and 100 fim flux densities of each resolution 
element then dividing them. We then sum the resolved luminosi- 
ties in our 4 different bands (thus L3_24, L24-70, L70-16O or L160-1100) 
and divide these integrated luminosities in each band by the inte- 
grated LjiR S500 (Table [TJ. This enables us to obtain the integrated 
fractions (y-axis). 

The thermal dust emission peaks in the [70-160 yt/m] band for 
all our galaxies. We thus naturally observe that the bulk of IR en- 
ergy is contained in that wavelength range. This [70-160 yt/m] frac- 
tion is homogeneous across the sample and accounts for ~40% of 
the hriR. The [24-70 fim] fraction ranges from 0.1 to 0.4, increas- 
ing with decreasing metallicity (or increasing 70/100 colour). Since 
low-metallicity environments usually contain warmer dust ( |Hunter| 
let al.|1989[rDale et al.|2005 ), the [24-70 jim] fraction is indeed ex- 
pected to significantly contribute to the total IR emission in those 
objects. 

While the [3-24 /um] band contributes to ~20% on average 
to the hriR for the sample, there is a large scatter in this frac- 
tion. This is likely associated with the contribution of Polycyclic 
Aromatic Hydrocarbons (PAH) emission. Several factors drive the 
scatter in the contribution of PAHs to the [3-24 //m] band. Low- 
metallicity galaxies usually show weak emission from PAHs ( |En-| 
Igelbracht et~ari [20051 [Jackson et aD '2006) and PAHs are also 
known to be sensitive to the hardness of the radiation field ( [Mad-[ 
|den et al.|20'06l|Gamano et al.|2003| [20051 ). The PAH size distri- 
bution also tends to be diff'erent at low metallicity, but whether the 
small PAHs are destroyed in the harsh conditions or whether they 
dominate because of diff'erent formation processes is still not clear 
Punt et al.[[20T0l [Sandstrom eTal. 2012J . The paucity of PAHs 
in low-metallicity environments could also be due to a delayed in- 
jection of carbon dust by AGB stars ( |Dwek| 19981 [Galliano et al.[ 
|2008| ). This low PAH emission is responsible for the weak [3-24 
yL/m] fraction in our metal-poor galaxies. Some metal-rich galaxies 
(NGC 5866, NGC 4594, NGC 1316) also present low [3-24yL/m] 
resolved fractions (resolved elements with L3_2a/^tir <15% and 
fy(70 yt/m)/fy(100 //m)<0.5). Those peculiar galaxies show only lit- 
tle dust emission relative to their stellar emission, as commented by 
IDraine et al.1 ( [2007l ). 

We observe a very clear trend of the submillimeter [160-1100 
fim] fraction with the 70/100 colour or the metallicity. The [160- 
1100 yL/m] band accounts for a few % for low-metallicity galaxies 
to up to 25% of the total infrared luminosity budget for metal-rich 
environments, consistent with the previous studies of [Dale et al.[ 
(|2001) and D ale & Helou[ ( [2002| . A log-lin fit of our data leads 
to the relation: Li6o-iioo/Lr/i? = -0.28 log(fv(70yL/m)/fv(100yum)) + 
0.07. 

While not shown, we also quantified the [500-1100 yt/m] 
fraction. This band contributes no more than 0.6% in the resolved 
elements of our sample, and only up to ~0.2% for the most 
metal-poor galaxies. Thus, even if present, we do not expect a 
submm excess to significantly modify our conclusions in this work. 



5 SPITZER AND HERSCHEL BANDS AS Ltir 
MONOCHROMATIC CALIBRATORS 

In this section, we want to study the resolved relationships linking 
the individual Spitzer and Herschel bands with the Ltir. A calibra- 



tion using the resolved elements of the complete sample of KING- 
FISH galaxies has been derived. However, we also analyse the re- 
lations for individual galaxies in order to get a handle on the scatter 
driven by the variety of our sample and study how individual rela- 
tions change with global galaxy characteristics or local ISM condi- 
tions (metallicity, 70/100 colour for instance). This will help us to 
unify the picture initiated by [Boquien et al.[ ( [201 1\ for M33 using a 
wider sample of galaxies, more representative of the diversity of the 
local Universe. We thus first analyse the individual and global rela- 
tions between monochromatic far-IR surface brightnesses and TIR 
surface brightnesses in Sections 5.1 and 5.2. We also derive empir- 
ical calibrations of Ljir from monochromatic far-IR luminosities 
in Section 5.3. We analyse the validity of our monochromatic cal- 
ibrations for the KINGFISH sample in Section 5.4. Finally, we in- 
vestigate in Section 5.5 the goodness of our calibrations for a wider 
range of environments, including high-redshift galaxies. 

5.1 A qualitative view of the relations 

We display the resolved TIR surface brightnesses (Stir) of our 
galaxies as a function of the IRAC 8 yum, MIPS 24 jim, PACS 70, 
100 and 160 jim and SPIRE 250 jim brightnesses (in W kpc'^, log 
scale, 3-cr detection) in Fig.|5] To study how metallicity influences 
those monochromatic relations, galaxies are colour-coded by 
metallicity (expressed as 12+log(0/H)), from metal-poor in red 
to metal-rich in dark purple colour. We overlay integrated values 
for comparison. Since we restrict our Ljir maps to pixels with 
a sufficient signal-to-noise in each Herschel bands, we could 
be missing flux in the faint outskirts of our objects. We prefer 
to use the global Spitzer and Herschel fluxes (from [Dale et aL] 
[2007 1 1 201 2 ) instead. These integrated fluxes are divided by the 
area covered by our selected pixels to obtain the average global 
monochromatic brightnesses in W kpc~^. In order to remove the 
stellar contribution from the integrated 8 jum flux densities of Dale| 
[et al.[ |2007 ), we apply the recipe of Marble etal.| 12010), namely 
fg^^'^^7f3.6~24%. This estimate is close to the value derived in 
[Helou et al.| ([2004) using the recipe of Starburst 99 (23.2%). We 
consider that the stellar contribution to the integrated 24 //m flux 
densities is negligible. 

We observe that : 

- the relation between the IRAC 8 jum surface brightness and 
Stir changes from galaxy to galaxy, as already shown by |Calzetti| 
[et al.| ( [2007| ), and strongly depends on the metallicity of the galaxy. 
For a given Stir, the resolved 8 yt/m brightnesses are systematically 
lower in metal-poor objects. This trend is consistent, as discussed 
in Section 4, with the decrease of PAH emission often observed in 
low metallicity environments, 

- MIPS 24 yL/m is also a decent tracer of Stir, but a significant 
scatter can be observed from galaxy-to-galaxy, as already observed 
by [Calzetti etaLlpOOT) ), 

- the relation between the brightness in the 70 fim band and 
Stir also seems to be slightly dependent on metallicity. For a given 
Stir, we observe higher 70 yt/m brightnesses in low-metallicity 
galaxies. The mid (MIR) to far-infrared (FIR) part of the SEDs of 
dwarf galaxies are known to be elevated compared to normal spi- 
ral galaxies, since they contain more small grains and hotter dust 
(Engelbrac ht et al.|2005[ [Galliano et al.|2005|or the globa l SEDs 
of KINGFISH low-metallicity objects in 'Dale et aL|20I2| ), which 
could explain this trend. The lowest metallicity still appears as an 
outHer in the PACS 160 and SPIRE 250 jim panel. 
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Figure 5. TIR surface brightness as a function of the individual bands (IRAC 8 //m, MIPS 24 jim, PACS 70 //m, PACS 100 //m, PACS 160 //m and SPIRE 250 
/im). Galaxies are colour-coded by metallicity (colour scale in the top left panel expressed as 12 + log(0/H)). We overlay global surface brightnesses with 
black circles. Dashed lines indicate the 1:1 relation. 



- the Herschel bands (PACS 70, 100 and 160 jdm bands, as 
well as SPIRE 250 jim band to a lesser extent) appear to be very 
good monochromatic indicators of ^tir as suggested by the tight 
correlations obtained in Fig. [5] 



We show the resolved Sjir (in W kpc~^, log scale) as a 
function of the different monochromatic surface brightnesses (in 
W kpc~^, log scale) for each galaxy of the KINGFISH sample in 
the Appendix (FigurejBTJ. We restrict this calibration to pixels with 
a 3-cr detection in the individual bands. We remind the reader that 
stellar contribution to the 8 fim and 24 fim bands is removed us- 
ing the stellar emission modelled during the SED fitting process. 
As mentioned before, this contribution is minor at 24 jim, with a 
median contribution of 0.67% over the pixels fulfilling our 3-cr cri- 
terion. 



5.2 Quantitative analysis 

For each galaxy and for the complete sample, we derive the cali- 
bration coefficients (a^, hi) such as: 



log Stir = log S i + bi 



(2) 



where Sjir refers to the TIR brightness, S/ the brightness in a given 
Spitzer or Herschel band / (from IRAC 8 yt/m to SPIRE 250 yt/m), 
both in W kpc~^, and a/ and b, respectively the slope and the inter- 
cept of the fit. 

Here and for the rest of Section 5, we choose to work in 
log-log space to account for non-linearities between the Ltir and 
the monochromatic emissions. We refer to |Boquien et al.| ( [201 1| | for 
a description of non-linearity eff'ects in the relations between TIR 
surface brightnesses and monochromatic surface brightnesses. We 
perform the regressions using a least-squares-bisector algorithm 
(function sixlin of the Astrolibrary of IDL) and use a "jack-knife" 
technique to quantify the goodness-of-fit and provide calibration 
coefficients with conservative errors. Indeed, we randomly select 
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N=l/10 of the resolved elements used for a given calibratiorj^ 
perform the regression for that subsample, save the coefficients 
and repeat this procedure 10000 times. The calibration coefficients 
provided in this paper are therefore the median of the coefficient 
distributions and errors on the coefficients are obtained from the 
standard deviations. 

We indicate the calibration coefficients obtained for the whole 
sample in Table |2] We also tabulate the coefficients a^ and b/ ob- 
tained for each galaxy in the Appendix in Table |C1| We add the 
scatter around the best fit (in dex) in the last column of the table. In 
order to better assess the results for individual galaxies, we plot in 
Fig.|6]the individual coefficients a/ and b/ derived using Eq.|2](and 
tabulated in Table [CTJ as a function of metallicity (expressed as 
12+log(0/H)). For galaxies with average 12+log(0/H) > 8.0, the 
MIR 8 and 24 fim brightnesses are linear estimators of the ^tir, 
as also found in Zhu et al. (2008 ). The three PACS bands are very 
reliable estimators of the Stir. The slope of the relation is, on av- 
erage, < 1.0, -1.0 and >1.0 for PACS 70, 100 and PACS 160 yt/m 
respectively, consistent with the slope we derive gathering the re- 
solved elements of the whole sample and tabulated in Table [2|. We 
observe that the scatter for the MIPS 24 fim and PACS relations is 
small, with a minimum scatter for the PACS 100 yt/m relation. 

For our low-metallicity galaxies, we observe that the 8 fim 
emission tends to be a sub-linear estimator of the Stir. This result is 
consistent with the trend observed for the resolved elements of low- 
metallicity galaxies in Fig. [5] Low-metallicity objects also seem to 
show steeper relations between PACS 100 and 160 jim monochro- 
matic brightnesses and Stir. This means that for a fixed PACS 
brightness, the TIR emission of low-metallicity environments will 
be higher than a normal spiral galaxy. We nevertheless highlight 
the difficulty in performing a pixel-by-pixel calibration from PACS 
data in some of our metal-poor galaxies due to poor statistics and 
a possible lack of detection at low-surface brightnesses that could 
bias the calibrations toward steeper relations. 

We finally note that residuals from the best fit to Eq.[2]for each 
galaxy are smaller than the residuals from the same fit performed 
on the integrated values of our galaxies. This could favour global 
parameters rather than local parameters as a driver of scatter in the 
relationship between monochromatic and TIR surface brightnesses. 

Comparison with M33 - [Boquien et al.| ( |2011| ) derived cal- 
ibration coefficients linking the monochromatic Spitzer and Her- 
schel surface brightnesses to Stir for the Scd galaxy M33, us- 
ing Eq. [2] We overlay the monochromatic coefficients they de- 
rive with red crosses in Fig. [6] for comparisons with our results 
for the KINGFISH galaxies. We adopt an oxygen abundance of 
12+log(0/H)=8.4 (from |Massey"[T998] ) for M33. For galaxies with 
12+log(0/H) > 8.0, the calibration coefficients derived in our study 
and in Boquien et a L| ([T0 1 1 ) for M33 are very similar. We note that 
[Boquien et aL| ( |2011| ) found that for the galaxy M33, the MIR 8 
and 24 fim brightnesses seem to be sub-linear estimators of Stir. 
To check if the difference could be linked with the signal-to-noise 
threshold we choose, we restrict our study to resolved elements 
above a higher (5cr) brightness threshold. This does not strongly 
modify our calibration coefficients for the 8 and 24 yum relations. 
A superimposition of the resolved elements of M33 with galaxies 



^ The choice of N does not influence our final results as long as the sub- 
sample still contains enough points to be representative of the relation. Co- 
efficients start to diff'er by a few percent if we choose N< 1/100 



of our sample sharing similar properties (metallicity, global Ltir, 
SFR) as M33 confirm the difference. Small differences could also 
arise as a result of differences in the regression methods. We thus 
apply the same linear regression (bisector algorithm) to the M33 
data used by [Boquien et al.| ( |2011| ) and obtain calibration coeffi- 
cients a8=0.88, b8=4.69, a24=0.84 and b24=6.61, very similar to the 
values found in their paper. The calibration coefficients of M33 are 
within the ranges of values found for our sample. Small differences 
could possibly due to 7 j a different treatment in the data reduction 
/ background subtraction at those wavelengths, 2j a real diff'erence 
in the 8 or 24 yt/m vs TIR surface brightness relations compared 
to KINGFISH galaxies with similar average oxygen abundance, 3) 
the large uncertainties on these average oxygen abundances, M33 
having a metallicity gradient for instance ( ,Magrini et al.|2007] ). 

5.3 Calibration from monochromatic luminosities 

Integrated luminosities are sometimes the only measure we can ac- 
cess, in particular for high-redshift observations. Studies of non re- 
solved objects thus require calibrations using luminosities rather 
than surface brightnesses. Gathering the resolved elements of the 
complete sample, we derive calibrations coefficients similar to 
those of Eq. [2]but linking Ltir and 24, 70, 100, 160 or 250 jum 
monochromatic luminosities. Eq.[2]can thus be re- written as: 

log Ltir = cii log vLy{i) + bi (3) 

with now Ltir the TIR luminosity, yLv(i) the flux in a given Spitzer 
or Herschel band / and both Ltir and yLy(i) are in L©. As before, we 
use our "jack-knife" technique to quantify errors on the calibration 
coefficients (see Section 5.2 for details on the technique) and report 
the coefficients calibrating Ltir from monochromatic luminosities 
in Table |2] The scatter around the best fit (in dex) is also provided 
in the table. 

If we were to derive a similar calibration for individual 
galaxies, the individual slopes (a^) would not be modified by the 
change of units compared to those derived in Section 5.2. Most 
of the conclusions of Section 5.2 still apply for the calibrations 
derived in this section. We note that the coefficients we derive 
for the 24 jim calibration (a24=0.954, b24= 1.336) are consistent 
with those derived by |Rieke et aL| ( [2009| ) within their error bars 
(namely a24=0.920, b24= 1.183). They estimate Ltir from IRAS 
fluxes using the recipe of |Sanders et"aLl ( [20031 ). 



5.4 Validity of the calibrations for KINGFISH galaxies 

We aim to test the robustness of the empirical calibrations derived 
in Section 5.3 - namely their ability to reproduce resolved and inte- 
grated luminosities - and analyse the intrinsic biases of our predic- 
tions. We thus compare the resolved and integrated Ltir calibrated 
from the monochromatic luminosities at 70, 100, 160 or 250 fim 
|3and the coefficients tabulated in Table [2] versus the resolved and 
integrated Ltir derived with a proper SED modelling. Figure |7] il- 
lustrates the comparison as a function of the 70/100 colour. For 
integrated values, we use the global flux densities of [Dale et al.| 
( |2012,) to predict integrated Ltir and compare them to our reference 

^ Because the PACS wavebands are close to the peak of the SEDs, they 
are expected to provide more robust calibrations of the Ltir than the 8 and 
24 //m wavebands. The 8 //m calibration is also strongly dependent on the 
metallicity. 
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Figure 6. Calibration coefficients to convert monochromatic brightnesses from Spitzer or Herschel bands into Stir plotted as a function of metallicity expressed 
as 12+log(0/H). The relation is : log Stir= a/ log S/ + b/ with S/ in W kpc~^. Horizontal dashed lines indicate a/=l (linear relation) and b/=0. We overlay the 
coefficients derived by ,Boquien et al.|{201H for M33 with red crosses and join an histogram of the coefficients on the right hand side of each plot. 



5500 (Table. [TJ. Results are identical for integrated or resolved 
hriR. The shift between resolved and integrated values (grey points 
/ filled black circles) may be linked with the fact that ( i) the cali- 
bration tree used in the data reduction of Herschel observations has 
changed since the Dale et al. p012j study and (ii) the background 
subtraction technique used in both studies are slightly diff'erent. 

PACS 70 - The 70 yt/m band provides a reasonably good 
monochromatic estimate of the Ltir (diff'erence < 50% for most of 



them). The Ltir of the lower IR luminosity objects (Ltir < 3 x 10^ 
Lq) is usually overestimated by the 70 yt/m monochromatic calibra- 
tion. Those objects are mostly low-metallicity galaxies (DDO 053, 
MSlDwB, Holmbergll, IC 2574, NGC 2915) that usually present 
warmer temperatures, which is consistent with the overestimation 
we observe. We also observe a strong correlation of the goodness 
of the 70 jim calibration with the 70/100 colour. 
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Table 2. Calibration coefficients to predict the TIR brightness/luminosity from monochromatic Spitzer and Herschel bright- 
nesses/luminosities. 



Surface brightnesses 
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Figure 7. Comparison between resolved and integrated I^tir calibrated from Herschel monochromatic fluxes I^tir band i (with i = 70, 100, 160 or 250 jim) 
with modelled I^tir 5500 plotted as a function of the fy(70//m)/fv(100//m) far-IR colour. The calibration relation is of the form log I^tir band i- a? log yLv(i) 
+ hi with luminosities in L©. We report the calibration coefficients in Table |Cl] (last line). We show the integrated values with black circles and the resolved 
elements with grey points. For each galaxy, we sort and average resolved element 10-by-lO for clarity. 
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PACS 100 - The 100 yt/m band offers the best monochromatic 
estimator, with very Uttle scatter and 53 out of 55 galaxies within 
~30% (as shown in NGC 6946 by Tabatabaei et al, submitted to 
A&A). We remind the reader that uncertainties on the resolved V.jir 
are of the order of 10-15%. We thus observe that many galaxies 
have predictions that match the modelled L^/^ within these error 
bars. The main outlier is NGC 1377. NGC 1377 shows a signifi- 
cant excess in its infrared-to-radio ratio and is thought to be a rare 
local nascent starburst probably powered by accretion through a re- 
cent merger ( [Roussel et al.|2006| . The galaxy presents a very hot 
IR SED peaking around 60 \im (yf^ units) and the R\CS 100 \im 
observation thus already belongs to the Rayleigh- Jeans slope of the 
SED of this object ( |Dale et a l. 2012). This hot SED compared to 
those of the other galaxies of the sample explains why we under- 
estimate its \.TiR when R\CS 100 \im is used as a monochromatic 
indicator (inability to capture the starburst component). The same 
explanation applies for the star-bursting dwarf irregular NGC 5408 
whose IR emission peaks at ~70 \m\ (Dale et al. 2012). Given the 
fact that the 70 \m\ band corresponds to the peak of the global SED 
for NGC 1377, using the 70 //m flux density combined with the co- 
efficients of the 100 \m\ calibration lead to a much better estimate 

of the \^TiR (hriR calibrated I ^TIR modelled - 0.76). 

PACS 160 - Like the 70 jjm calibration, the goodness of our 
160 iim calibration depends on the 70/100 colour. The hot SED of 
NGC 1377 or NGC 5408 can explain why our 160 yt/m monochro- 
matic calibration underestimates the integrated I^tir for these ob- 
jects. DDO 053 is also an outlier for the R^CS 160 iim calibration. 
DDO 053 shows a very low integrated PACS 160 imv flux compared 
to the PACS 100 and SPIRE 250 fim fluxes (see'Dale et al.»20121). 
Using the MIPS 160 jim flux of this galaxy estimated in Dale et al. 
([2007} (0.5 Jy in Ueu of 0.25 Jy for PACS 160 //m) would lead to 
a better agreement of the calibrated TIR value with the modelled 
integrated I^tir- 

SPIRE 250 - We finally observe a larger uncertainty in the 
predictions derived using the 250 yum monochromatic calibration 
(still within a factor of 2 for most of our objects) with, as expected, 
a strong correlation of the goodness of the calibration with the 
70/100 colour similar to that of the 160 yum calibration. Our 
calibration will thus probably underestimate the Ljir for hot 
objects and overestimate the Ltir for cold objects. 

We conclude that PACS 100 yum luminosities can be safely 
used as monochromatic estimators of the Ltir, even if it should be 
used with caution for strong starburst environments. We note the 
small errors on the parameters derived and small scatter around the 
relation. Our 70 and 160 /um calibrations also lead to reasonably 
good estimates of the Ltir (within 50% for most resolved elements 
or entire galaxies). Here again we caution their use for very cold 
or very hot SEDs. Calibrations using the 250 fim luminosity alone 
have larger uncertainties. We observe that predictions deviate from 
modelled values if fvo^fioo or if f7o/fioo<0-4. Calibrations using 
combined fluxes should be favored if more than one PACS/SPIRE 
band is available (see Section 6). 

In the present paper, we estimate calibrations of the Ltir from 
observations in the [8-250 fim] band, but calibrations at longer 
wavebands would be useful to understand galaxy properties of 
nearby or high-redshift objects observed in submm and millimeter, 
from ground-based telescopes in particular (SCUBA-2, LAB OCA, 
ALMA). In Appendix A, we present and discuss monochromatic 
calibrations derived for longer wavelengths from 1 ) 350 and 500 
Herschel data 2) 850 and 1000 jim model predictions. 



5.5 Predictions for near and high-redshift sources 

We now analyse Ltir predictions for a wider range of environ- 
ments using our monochromatic calibrations. Figure [8] compares 
the Ltir modelled in this paper or published in the literature of 
various objects with their Ltir predicted from our monochromatic 
calibrations (and their rest-frame 70, 100, 160 and 250 yt/m lumi- 
nosities). Black points indicate the resolved elements of the KING- 
FISH galaxies, with integrated values overlaid with black circles. 
The horizontal dashed line indicates when predictions from our cal- 
ibrations match modelled or published Ltir. We add various nearby 
and high-redshift objects for comparison. 



5. 5. 1 Nearby sources 

|Le Floc'h et al.| 12012) characterise the close environment of the 
gamma-ray burst GRB 980425, located 36 Mpc away. They pro- 
vide 70 and 160 yt/m fluxes for the GRB host (respectively 230 and 
615 mJy). Using our monochromatic calibrations, we predict an in- 
tegrated Ltir of 8.17^0-6 x 10^ and l.OT^-^ x 10^ from the 70 
and 160 fim fluxes respectively. Those estimates are in very good 
agreement with the integrated Ltir they derive using standard em- 
pirical libraries of galaxy templates (1.02 x 10^ L©). We overlay 
the environment of GRB 980425 in Fig. [8] (70 and 160 yt/m panels, 
upside-down orange triangles). 

|Rangwala et al.| |201 1) present SPIRE-FTS Herschel observa- 
tions of the nearby ultra-luminous infrared galaxy Arp 220, located 
at 77 Mpc, and provide a SPIRE 250 jim continuum flux of 30.1 
Jy for this galaxy. Using a single-temperature modified blackbody 
to fit the global SED from 15 fim up to SPIRE bands (including 
ISOPHOT, IRAS, ISO-LWS, SPIRE and SCUBA data), they esti- 
mate a Ltir of 1.77 x 10^^ L©. From our 250 yum monochromatic 
calibration, we estimate the integrated Ltir to be 7.95-^-^ x 10^^ L©, 
thus a factor of two lower than their modelled value (see Figure [8] 
250 jim panel, upside-down purple triangle). Our 250 jjm calibra- 
tion is expected to under-predict the integrated Ltir for hot objects 
(Arp 220 has an IRAS 60-to-lOO ratio ~ iQ which could explain 
the diff'erence. Assuming the 100 yum flux is 130 Jy (from ISO- 
LWS), we derive a Ltir of 1.23^^-^ X 10^^ f^j. ^^is object, thus 
~70% of value derived by Rang wala et al . | pO^ 11) . This highlights 
the good predictions from our 100 yum monochromatic calibration, 
especially for unusual object such as Arp 220. 

Finally, Pereira-Santaella et al. (in prep) present SPIRE-FTS 
Herschel observations of local active galaxies and provide SPIRE 
250 yL/m continuum flux densities for those sources. They also pro- 
vide integrated IR luminosities taken from [Sanders et al.| p003| ) 
and rescaled to the distance they adopt. For most of the galaxies of 
their sample, the SPIRE fluxes only include the nuclear far-IR emis- 
sion while the IR luminosities are derived from the integrated IRAS 
fluxes. Comparisons should be safe for UGC05101 and NGC 7130. 
We consider the other 250 fim fluxes as lower limits. We overplot 
these objects with red upward arrows in Figure |8] (250 jjm panel) 
and UGC05101 and NGC 7130 with red crosses. Our predictions 
using the 250 yt/m calibration match the IR luminosities published 
in their study for these two objects within a factor of 2. 



^ IRAS flux densities taken from the NED database 
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Figure 8. Comparison between 'Ltir (modelled in this paper or published in the literature) of various objects and their I^tir predicted from our monochromatic 
calibrations (and their rest-frame 70, 100, 160 and 250 iim luminosities). Black points indicate the resolved elements of the KINGFISH galaxies. For each 
galaxy, we sort and average resolved elements 20-by-20 for clarity. We overplot integrated measures of our KINGFISH sample with black circles. Orange 
upside down triangles indicate the environment of the gamma-ray burst ORB 980425 from jLe Floc'h et aT]j2012) . Magenta triangles indicate submillimeter 
galaxies from Magn elli et ar]j2012) . We also overlay the uncertainties on their modelled I^tir- Green circles indicate z~2 dust obscured galaxies from 
[Melbourne et al. (20 12). Yellow diamon ds indicate WISE-Selected Hyper-luminous Galaxies from |Wu et al.H20l'2) . The blue star shows an Hyper-Luminous 
Infrared Galaxy from |Eisenhardt et al.H201 2). Red crosses and upward arrows (lower limits) indicate local active galaxies from Pereira-Santaella et al (in 
prep). The upside down purple triangle finally shows Arp 220 studied ,Rangwala et aT]j2011) . 



5. 5. 2 High redshift sources 

We combine these nearby results with a selection of z~l-3 objects 
taken from published catalogs. We derive distances from the red- 
shifts provided in the catalogs, adopting an = 70 km s~^Mpc~\ 
Qm = 0.3, Qa = 0.7 cosmology. In each catalog, we select objects 
observed at rest-frames close to 70, 100, 160 or 250 jim, within a 
margin of +10 fim, and compare the published integrated I^tir with 
those derived from our calibrations. For each sample, we quote the 
methods used to derive the published I^tir to help the reader assess 
the level of confidence of these modelled values and try to under- 
stand how much of the scatter in Fig.[8]could be driven by biases in 
our predictions or, on the contrary, linked with uncertainties in the 
published values we quote (derived with dififerent methods and, for 
some of them, with limited data). 

We first select z~2 dust obscured galaxies from the catalog 
of [Melbourne et al.] ( |2012| l (Fig. [8] green circles). Those galaxies 
are observed with Herschel at 250, 350, 500 iim. They estimate 
the integrated I^tir by interpolating between the mid to far-IR flux 
densities and extrapolate the long wavelength tail contribution to 
the l^TiR using a modified blackbody curve, assuming a dust emis- 



sivity index of 1.5. Even though they do not completely sample the 
submm slope of their high-redshift sources, they note that their es- 
timated l^TiR are relatively robust, with less than a 5% change in 
l^TiR if the temperature of the modified blackbody they used varies 
by 25%. In Fig.jS] we observe that our 70, 100 and 160 /im pre- 
dictions are consistent with the integrated I^tir estimated by |Mel-[ 
[bourne et al.[ ( [2012t within a factor of 3. We remind the reader 
that we select objects observed at rest-frames close but not exactly 
equal to 70, 100, 160 or 250 jim. This results in some cases in an 
under or over-prediction of the monochromatic flux used in the cal- 
ibration and thus of the predicted I^tir- Restricting our selection to 
a margin of +5% around our reference wavebands remove some 
of the objects showing high discrepancies. This does not, however, 
fully explain the vertical scatter we obtain. The galaxies for which 
our 100 jim calibration under-predicts the I^tir are all classified as 
Mrk231-like objects (namely AGN-dominated ULIRGs) in |Mel-[ 
bourne et al. (2012), with temperature superior to 40K. The dis- 
crepancy between the modelled and our predicted I^tir can thus be 
attributed to the fact that these objects have SEDs that peak at much 
shorter wavelength than 100 jim (Section 5.4). Stronger discrepan- 
cies are observed when using the 250 iim calibration. For the three 
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z~l galaxies detected at 500 jim (so rest-frame 250 jim), namely 
J143052.8+342933, J143313.4+333510 and J143334.0+342518, 
our 250 jim calibration predicts a I^tir higher by a factor of 5.5 
to 8.7 compared to those derived in Melbour ne et^L] ( [2012| ). Those 
galaxies are among the coldest objects of the sample (T<26K). Our 
250 jim calibration probably over-predicts the I^tir in those objects. 

We also add submillimeter galaxies (SMGs) from |Magnelli| 
letaTlpOTIl ) observed with Herschel at 160, 250, 350 and 500 (im 
(magenta triangles). The published integrated I^tir are derived us- 
ing a power-law temperature distribution model. We include the 
lensed-SMGs in our analysis. Herschel flux densities of these ob- 
jects have been de-magnified using the magnification factors tab- 
ulated in the Table 11 of |Magnelli et al.| { [2012| ). Uncertainties on 
the modelled I^tir are provided for this sample and were added 
to Fig. [8] We observe that the I^tir of submillimeter galaxies esti- 
mated in [Magnelli et al.| ( |2012| ) are systematically lower than our 
predictions, whatever the rest-frame used. Our monochromatic 70, 
100 and 160 fim calibrations predict a I^tir higher by a factor of 3 
at most compared to those derived by iMagnelli et al.| ( |2012| ), our 
250 jim calibration by up to an order of magnitude. For the top 
left panel (70yL/m), deviations from our predictions seem to increase 
with luminosities. As reminded in "Magnelli et al.| ( [2012| ), the SMG 
population is very heterogenous and biased towards cold dust tem- 
peratures compared to the entire infrared galaxy population. This 
can partly explain why the I^tir predictions using our 160 or 250 
jjm calibrations are higher than the modelled values. If the SED 
profiles of SMGs were similar in the shape than those of local ob- 
jects but simply shifted to shorter wavelengths, the inverse trend 
would be observed at 70 jim, namely that our calibration would 
under-predict the I^tir- This is not what is observed in Fig.jS] We 
finally note that part of the TIR luminosity could be missing in the 
modelled I^tir values of Magnelli et al.'('2012 ) due to the fact that 
their fits do not include a detailed modelling of the rest-frame [8 to 
70 iim\ spectrum. Indeed, SMGs often show broad emission fea- 
tures from PAHs and starburst activity could dominate the I^tir in 
those objects (e.g. |Menendez-Delmestre et al.|2009| ). 

We finally add WISE- selected hyper-luminous galaxies from 
|Wu et al.| ^Yl) observed with CSO/SHARC-II at 350 and 
450 fim and CSO/Bolocam at 1.1 mm (yellow diamonds) and 
a z=2.452 hyper-luminous galaxy from |Eisenhardt et al.| ( |2012| ) 
(WISE 1814+3412) observed with CSO/SHARC-II at 350 iim 
(blue star, 100 fim panel). The same SED fitting technique is ap- 
plied in the two studies to derive integrated TIR luminosities: they 
use a single modified blackbody model combined with power-laws 
to connect the mid-IR to mm SED points and a modified blackbody 
component (with a dust emissivity index of 1.5) to fit the longer 
wavelengths. Their integrated I^tir can thus be considered as low 
limits for the total luminosities. The modelled I^tir values from 
|Wu et al.| ( [20T2l l and |Eisenhardt et al.| { [20T2j ) match with our 100 
^m calibration, here again within a factor of 3. |Wu et ar] ( |2012| ) 
note that their objects could host highly obscured AGNs heating 
their dust cocoon to very high temperatures. This could explain 
why our 100 jim calibration under-predict the I^tir- We also note 
that in the fourth pa nel of Fig. [8) our predicted I^tir for the galaxy 
W0149+2350 (from |Wu et al.|2012| ) is lower than their modelled 
value by a factor of 7.6 while our prediction using the 100 jdm cal- 
ibration match their modelled value within 30%. This galaxy has a 
redshift of z=3.228 and was observed with the Submillimeter Ar- 
ray (SMA) at 1.3 mm. The 1.1mm flux density quoted in |Wu et aT] 
( |2012| ) is a conversion of the SMA measurement assuming an emis- 
sivity index and 2, then taking the average of the two values. 



Discrepancy could thus be linked with the uncertainty on the pre- 
dicted 1.1mm flux for this object. 

We conclude that while we observe a wide spread for objects 
with integrated I^tir > 10^^ L© on the 250 jim panel, predictions of 
the PACS monochromatic calibrations match modelled I^tir from 
resolved elements of nearby galaxies to global objects at further 
redshift over a surprisingly large luminosity range (from 10^ to 
10^^ Lq). This reinforces the usefulness of PACS wavelengths as 
reliable monochromatic calibrators for the I^tir of nearby galaxies. 
This also means that reliable estimates of the I^tir of high-redshift 
objects can be obtained using the SPIRE filters and our 70, 100 
and 160 jjm calibrations. For instance, SPIRE 250 jim observations 
combined with our 100 jim calibration coefficients could be used 
to study the peak of the star formation history at z=1.5. 



6 COMBINING MID TO FAR-IR BANDS 

[Dale & Helou| ( |20()2] ) derived a bolometric relation to estimate the 
integrated I^tir from a combination of MIPS filters. This relation 
is applicable to a wide range of galaxy luminosities. With Her- 
schel observations, our wavelength coverage now goes longward 
of 160 fim. This enables us to better sample the submm slope of the 
SED and reduce the uncertainties on integrated I^tir linked with 
the presence of cold dust not detected by previous MIPS 160 jim 
observations. In this section, we thus aim to make the most of the 
good resolution of Herschel to perform a similar multi- wavelength 
empirical calibration on a resolved basis. We derive calibrations 
linking combined far-IR brightnesses (Section 6.1) or luminosities 
(Section 6.2) to the TIR surface brightnesses or luminosities and 
study their dependence and biases. 



6.1 Calibration from combined brightnesses 

Using a combination of resolved brightnesses at MIPS 24 yt/m, 
PACS 70, 100 and 160 fim and SPIRE 250 yum, we calculate caU- 
bration coefficients C/ for each galaxy, such as: 

^r/i? = ^Q5, (4) 

where Stir refers to the TIR surface brightness and S^ the bright- 
ness in a given Spitzer or Herschel band /. Here and for the rest of 
Section 6, we derive calibrations in linear space. We use a multi- 
ple linear regression fit (function mregress, a variant from the IDL 
function regress by Ph. Prugniel, 2008) combined with our "jack- 
knife" technique applied on the resolved elements of our galaxies to 
conservatively estimate the calibration coefficients and their uncer- 
tainties. We list the coefficients obtained for individual galaxies in 
the Appendix in Table |C2| We indicate the calibration coefficients 
obtained for the whole gathering the resolved ISM elements of the 
whole sample in Table [3] 

To quantify the scatter between the modelled and the predicted 
brightnesses for each combination, we also provide indicators of 
the goodness-of-fit in the last two columns: the coefficient of de- 
termination R^ and the coefficient of variation of the root-mean- 
square error CV(RMSE). R^ ranges between and 1 and indicates 
the proportion of variability of the resolved TIR brightnesses ac- 
counted for by our calibration. For instance, R^=0.90 means that 
our calibration accounts for 90% of the total variation of our TIR 
brightnesses. CV(RMSE) is the standard deviation (measuring the 
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Table 3. Calibration coefficients to predict the TIR brightness/luminosity from combined Spitzer and Herschel bright- 
nesses/luminosities. 
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Figure 9. Calibration coefficients to derive Stir from a combination of 
Spitzer and Herschel brightnesses (24, 70, 100, 160 + 250 //m) plotted as a 
function of metallicity expressed as 12+log(0/H). The relation is : Stir = Z 
Ci Si with Stir and S/ in W kpc~^. We join an histogram of the coefficients 
on the right hand side of each plot. 



differences between the TIR brightnesses predicted by our calibra- 
tions and the TIR brightnesses we obtained using the [DL07] mod- 
elling) normalised to the mean values of our resolved TIR bright- 
nesses. Therefore, the lower the CV(RMSE), the better. These two 
quantities are defined as: 



1 



CV(RMSE) = 



E (Mi - Pif 
I.(Mi - Mif 

RMSE 



Mi 



1 

Mi 



^{Mi-Pif 



(5) 



(6) 



modelled brightnesses and n the number of ISM elements. 

We plot the coefficients C/ derived using Eq. |4] as a function 
of metallicity in Fig. [9] Calibration coefficients weighting the 
PACS 70 and 100 fim brightnesses are similar though the sample 
as suggested by the peaked distributions of the histograms. We 
nevertheless observe a larger distribution in the 160 and 250 fim 
coefficients. No strong trend is observed with metallicity for the 
24, 70, 100 and 160 jjm coefficients. However, the C250 coefficients 
are quite low for objects with 12+log(0/H) < 8.2 (median=-0.01), 
probably linked with he fact that submm emission has a smaller 
contribution to the total IR budget in low-metallicity galaxies 
compared to more metal-rich objects. 



6.2 Calibration from combined luminosities 

Similarly to Section 5.3, we gather the resolved elements of the 
complete sample and derive calibration coefficients C/ similar to 
Eq.[4]but linking the I^tir to different combinations of Spitzer and 
Herschel luminosities. Eq.[4]can thus be re- written as: 



^ Ci vLy{i) 



(7) 



with the predicted surface brightnesses, the resolved TIR 
surface brightnesses modelled using [DL07], Mi the mean of the 



where I^tir now refers to the TIR luminosity and yLy(i) the 
resolved luminosities in a given Spitzer or Herschel band /. I^tir 
and the different yLv(i) are in L©. We list the calibration coefficients 
obtained for various combinations of Spitzer and Herschel bands 
in Table[3]with respective R^ and CV(RMSE) coefficients. 



6.3 Validity of the calibrations for KINGFISH galaxies 

Following the same scheme as Section 5, we test the ability of our 
combined calibrations (derived using all the resolved elements of 
KINGFISH galaxies) to predict resolved and integrated I^tir and 
analyse their intrinsic biases. We thus compare the resolved and 
integrated I^tir of the KINGFISH sample calibrated from various 
combinations of the 24, 70, 100, 160 or 250 yum luminosities ver- 
sus the l^TiR derived with a proper SED modelling. For integrated 
values, we use the integrated flux densities of |Dale et al.|p007l 
|2012| | to predict integrated I^tir and compare them to the modelled 
integrated I^tirssqo (Table [T|. Figure [To| illustrates these compar- 
isons as a function of the 70/100 colour. Results are identical for 
integrated or resolved I^tir- 

We do observe a correlation between the goodness of some 
calibrations ( "24+70", "160+250", "100+160+250" for instance) 
and the 70/100 colour. Nevertheless, our combined calibrations 
lead in most cases to better estimates of the I^tir than monochro- 
matic calibrations. Combined calibrations are especially much 
more reliable than our 70, 160 or 250 pm monochromatic calibra- 
tions for galaxies whose 70-to-lOO flux density ratios are below 
0.4 or above 0.8. We note that combining MIPS 24, PACS 70 and 
PACS 160 pm data alone, we obtain C24 = 2.126, C70 = 0.670 and 
C160 = 1.134. We remind the reader that the coefficients obtained by 
[Dale & Hel6u| (p002 ) and calibrated from Spitzer 24, 70 and 160 
pm fluxes are 1.559, 0.7686 and 1.347 for C24, C70 and Ci6o A^m re- 
spectively. Global l^TiR derived using the Pale & Helou|(2002| ) cal- 
ibration differ from our modelled I^tir 5500 by ~26% at most (for 
IC 2574) with a median of the differences of ~6%. The integrated 
values obtained with the "24+70+160" pm calibration presented in 
this paper differ from the modelled integrated I^tir by ~ 23% at 
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Figure 10. Comparison between the resolved (grey points) and integrated (black circles) Ltir obtained with a combination of monochromatic fluxes (among 
24, 70, 100, 160 and 250 jum) and the properly modelled Ltir as a function of the fv(70//m)/fv(100//m) far-IR colour. The calibration relations are of the form 
Ltir= S Ci yLy(i) with Ltir and vLy{i) in L©. We indicate the 2, 3, 4 or 5 bands used to obtain the calibrated TIR luminosities in each panel. We report the 
various calibration coefficients in Table[3] We sort and average the resolved elements (in grey) 15-by-15 for clarity. Note that the Y-axis range is diff'erent from 
that of Fig.|7] 
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most (for NGC 1512), with a median of the differences of ~4%. 
Both predictions are thus very close at global scale (integrated lu- 
minosities). 

The calibration performed using the complete coverage of the 
FIR / submm emission (24, 70, 100, 160 and 250 fim data) leads 
to a reliable approximation of the modelled I^tir (low CV(RMSE) 
value) : using this calibration, all the estimated luminosities reside 
within 19% of the modelled I^tir, with a median of the differences 
of 3.5%. Using 4 of these wavelengths leads to similarly good re- 
sults. Any l^TiR predicted from a fewer number of fluxes should 
contain the 100 fim flux, or a combination of 70+160 to lead to 
l^TiR predictions reliable within 25%. We note that calibrations in- 
cluding 24 yL/m data lead to better estimates of the I^tir for galaxies 
showing high 70/100 colour than calibrations that do not include 
this wavelength. We thus advice using the 24 jim flux (if available) 
in the I^tir predictions for these environments. 



7 CONCLUSIONS 

• We investigate how SPIRE wavelengths influence the de- 
termination of the l^TiR and conclude that using data up to 250 
jjm leads to I^tir values that are in very good agreement with 
that obtained with a complete SED modelling of the dust thermal 
emission (within 10% for most of our resolved elements). 

• The [70-160] band contains 30 to 50% of the IR emission. 
We observe an overall shift in the SED to shorter wavelengths 
with decreasing metallicity. Indeed, the [24-70 fim] fraction 
increases for warmer sources (often found in low-metallicity 
objects) while the [160-1100 jim] fraction accounts for only a few 
% for low-metallicity galaxies (to up to 25% of the total infrared 
luminosity budget for metal-rich environments). The [3-24 //m] 
fraction accounts for ~20% of the L^/^, with a significant scatter 
from one environment to another. 

• We study the correlation between TIR and monochromatic 
Spitzer and Herschel surface brightnesses/luminosities and derive 
calibration coefficients to quantify these correlations. For most of 
the galaxies of our sample, the three PACS bands can be used as 
reliable monochromatic estimators of L^/^, with slopes on average 
< 1.0, -1.0 and > 1.0 for 70, 100 and 160 yt/m respectively. We also 
observe a strong correlation between the SPIRE 250 jim and L^/^, 
although with more scatter than the PACS relations. We estimate 
calibration coefficients for waveband beyond 250 fim in Appendix 
A. 

• We conclude that the 100 fim band is the best band to use as a 
monochromatic estimator (scatter of 0.05 dex) of L^/^. 

• We show that the calibrations at 70, 100 and 160 fim re- 
produce modelled I^tir over a very large luminosity range, from 
nearby galaxies to galaxies at z~l-3. I^tir values are reproduced 
with larger uncertainties from 250 yt/m fluxes. We nevertheless 
caution the use of our 70, 160 and 250 fim calibration for strong 
star-bursting environments and for objects showing cold dust 
temperatures. 

• We finally derive calibration coefficients to derive TIR surface 
brightnesses/luminosities from a combination of Spitzer and Her- 
schel surface brightnesses/fluxes. These calibrations lead to better 
estimates ofl^rm than monochromatic calibrations and show much 



smaller biases. We update the widely used I^tir calibration of |Dale| 
|& Helou| ( [2002t using Herschel/FACS 70 and 160 yt/m data in Heu 
of Spitzer/MIFS data at the same wavelengths. The two calibrations 
lead to similar estimates (with similar uncertainties) for integrated 
luminosities. As expected, the calibration using the complete sam- 
pling of the FIR/submm emission (data at 24, 70, 100, 160 and 
250 yL/m) leads to a reliable estimation of the Ltir but using 4 of 
those wavelengths leads to similarly satisfying predictions. We note 
that including 24 fim data in the calibration is essential to properly 
estimate Ltir in strongly star-forming environments (high 70/100 
colour). 
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at those wavelengths (free-free or synchrotron emission, molecu- 
lar line contamination etc) has to be removed if the calibration is 
used. The slopes of the 500, 850 and 1000 jjm calibrations are very 
close because these observations sample the submm slope of the 
SEDs where fluxes are evolving in a similar way. The 850 and 1000 
maps are moreover directly extrapolated from the SED model per- 
formed using data up to 500 yum. Our results are thus consistent 
with a scaling of the fluxes (Lgso/Liooo^constant) from one relation 
to the other, translated in log space by a simple shift of the intercept. 
Small differences are nevertheless due to the fact that the submm 
slope varies from one galaxy to another. 



APPENDIX A: TIR CALIBRATION BEYOND 250 yt/m 
Al 350 and 500 jim 

To estimate similar monochromatic calibrations for the 350 and 500 
jjm wavebands, we use the I^tir maps obtained at SPIRE 500 res- 
olution. The FWHM of the PSF for SPIRE 500 fim is -36'' and 
the pixel size of our I^tir maps is lA" . We refer to Section 3.1 for 
explanations of the methodology. We gather the resolved elements 
of the complete sample to derive calibrations coefficients linking 
the resolved TIR luminosities with the 350 or 500 //m fluxes. We 
remind Eq.|3]here: 

log Ltir = at log yLy(i) + bf (Al) 

where Ltir refers to the TIR luminosity and v Ly(i) the flux at 
350 or 500 fim. Both Ljir and v hy(i) are in L©. We refer to Section 
5 for details on the regression technique. 

We obtain the following calibration coefficients : 

• ( a35o , b35o ) = ( 1-106 + 0.011 , 0.661 + 0.064 ) 

• ( asoo , bsoo ) = ( 1-160 + 0.012 , 1.008 + 0.062 ) 

The 350 and 500 /um relations with the TIR luminosities are 
thus over-linear relations. As for the 250 jum band (Fig.|5]), we ob- 
serve a large spread of Ljir values for a given 350 or 500 yum flux 
as well as a strong correlation of the goodness of our 350 and 500 
fim monochromatic calibrations with the 70-to-lOO flux density ra- 
tio. We thus caution the use of these calibrations for extreme (star- 
forming or very cold) environments. 



APPENDIX B: TIR SURFACE BRIGHTNESS VERSUS 
SPITZER/HERSCHEL BRIGHTNESSES 



A2 850 and 1000 jim 

From our resolved SED modelling performed at the resolution 
of SPIRE 500 Jim, we also extrapolate maps of the KINGFISH 
galaxies at 850 and 1000 yum, wavebands that are frequently 
observed from ground-based telescopes (SCUBA-2, LAB OCA, 
MAMBO etc.). We gather the resolved elements of the complete 
sample and estimate calibration coefficients to link the resolved 
TIR luminosities to the extrapolated 850 or 1000 yt/m resolved 
fluxes, both in L© (Eq.[3]). 

We obtain the following calibration coefficients : 

• ( agso , bgso ) = ( 1-150 + 0.013 , 2.161 + 0.057 ) 

• ( aiooo , biooo ) = ( 1-152 ± 0.014 , 2.533 ± 0.054 ) 



We remind the reader that the extrapolated 850 and 1000 jim 
fluxes used in this calibration are coming from pure thermal dust 
emission. Any non-dust contamination contributing to observations 
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Figure Bl. TIR surface brightness as a function of the different Spitzer and Herschel bands for the galaxies of the KINGFISH sample. We overlay the 
regressions estimated for each band. For each galaxy, we sort and average resolved element 5-by-5 for clarity. We report the coefficients of the fits in Table |Cl| 
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Figure Bl. continued. 
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Table CI. Calibration coefficients derived for individual galaxies to convert the monochromatic brightnesses in various Spitzer 
or Herschel bands (8, 24, 70, 100, 160 or 250 iim) into TIR brightness. 



Name 


ag bg 




a24 


^24 




a70 






DDO053 





630± 


0.043 


13.339± 1.390 





771± 


0.028 


8.350+ 0.900 


1 


031± 


0.219 


-0.779± 7 


.364 


Holmbergll 





718± 


0.028 


10.492± 0.895 





623± 


0.011 


13.238+ 0.352 


1 


161± 


0.019 


-5.243+ 


.638 


IC342 


1 


079 ± 


0.004 


-2.024+ 0.131 





825± 


0.003 


6.926± 0.096 





916± 


0.003 


3.293+ 


106 


IC2574 





613± 


0.031 


13.926± 0.999 





653± 


0.051 


12.433± 1.661 





899± 


0.053 


3.542+ 1 


801 


M81Dw 


1 


210± 


0.120 


-5.409± 3.827 





705 ± 


0.094 


10.675± 3.050 





724± 


0.029 


9.385± 


966 


NGC0337 





945+ 


0.010 


2J5o± 0.343 





896+ 


0.008 


4.563+ 0.278 





998+ 


0.009 


0.334± 


301 


NGC0628 





992+ 


0.006 


0.917+ 0.205 





881+ 


0.007 


5.000+ 0.247 





904+ 


0.007 


3.672+ 


229 


NGC0855 





914+ 


0.012 


3.836+ 0.386 





881+ 


0.008 


5.069+ 0.258 





972+ 


0.006 


1.177+ 


198 


NGC0925 





812± 


0.008 


7.101± 0.279 





867± 


0.008 


5.517± 0.255 





943 ± 


0.006 


2.273± 


198 


NGC1097 


1 


103± 


0.005 


-2.749± 0.184 





904± 


0.003 


4.259± 0.111 





875± 


0.004 


4.702± 


136 


NGC1266 


1 


021± 


0.027 


0.746± 0.898 





977± 


0.023 


1.699± 0.781 





965± 


0.015 


1.329± 


526 


NGC1291 





949± 


0.020 


2.655+ 0.671 





817± 


0.017 


7.186± 0.557 





864± 


0.009 


4.895± 


318 


NGC1316 


1 


076± 


0.015 


-1.573± 0.483 





965± 


0.024 


2.318± 0.788 





849± 


0.009 


5.435± 


297 


NGC1377 





963± 


0.008 


1.938± 0.270 





957± 


0.005 


2.018± 0.183 


1 


008± 


0.005 


0.075± 


174 


NGC1482 


1 


036± 


0.011 


-0.516± 0.372 


1 


081+ 


0.012 


-1.776± 0.399 





961± 


0.004 


1.619± 


139 


NGC1512 





997+ 


0.019 


rv oo'~i 1 rv 

0.887+ 0.637 





821+ 


0.017 


7.018+ 0.565 





865+ 


0.008 


4.961+ 


284 


NGC2146 


1 


009+ 


0.006 


rv orv /I 1 rv 1 rvc 

0.394+ 0.195 


1 


019+ 


0.005 


0.318+ 0.181 





952+ 


0.003 


1.943+ 


116 


NGC2798 


1 


033+ 


0.013 


rv T70 1 rv Aoo 

-0.278+ 0.438 


1 


094+ 


0.018 


-2.271+ 0.609 





948+ 


0.005 


2.021+ 


188 


NGC2841 


1 


090± 


0.010 


-2.312+ 0.352 





963± 


0.008 


2.371+ 0.281 





793± 


0.023 


7.540± 


778 


NGC2915 





916± 


0.021 


3.883± 0.695 





819+ 


0.015 


7.058± 0.486 





93 1± 


0.009 


2.560+ 


297 


NGC2976 





993± 


0.008 


rv rvrvo . rv '~\ z' f\ 

0.998± 0.269 





895+ 


0.008 


4.627± 0.262 





916± 


0.005 


3.227± 


167 


NGC3049 


1 


036± 


0.023 


rv o crv ■ rv tio 

-0.350± 0.778 





861 + 


0.013 


5.508± 0.432 





966± 


0.012 


1.451± 


421 


NGC3077 


1 


059± 


0.009 


1 /I . rv o rvrv 

-1.224± 0.309 





934± 


0.007 


3.237± 0.242 





912± 


0.005 


3.314± 


175 


NGC3184 


1 


080± 


0.015 


f\c c 1 rv /I rv 1 

-2.055± 0.491 





835+ 


0.008 


6.564± 0.265 





855± 


0.007 


5.355± 


249 


NGC3190 





961+ 


0.018 


2.153+ 0.593 





974+ 


0.014 


2.109+ 0.447 





947+ 


0.021 


2.183+ 


729 


NGC3198 





854+ 


0.012 


5.604+ 0.387 





806+ 


0.007 


7.485+ 0.237 





975+ 


0.010 


1.288+ 


328 


NGC3265 


1 


008+ 


0.016 


0.520+ 0.528 





951+ 


0.010 


2.445+ 0.334 





998+ 


0.013 


0.342+ 


443 


NGC3351 


1 


215± 


0.014 


-6.456± 0.457 





853± 


0.006 


5.927± 0.197 





852± 


0.007 


5.425± 


254 


NGC3521 


1 


005+ 


0.002 


0.461 ± 0.076 





947+ 


0.002 


2.889± 0.072 





895+ 


0.002 


4.039± 


072 


NGC3621 





940± 


0.004 


2.687± 0.134 





939+ 


0.004 


3.152± 0.121 





916± 


0.004 


3.257± 


132 


NGC3627 


1 


028± 


0.005 


-0.236± 0.170 





984+ 


0.005 


1.618± 0.157 





885± 


0.004 


4.372± 


143 


NGC3773 


1 


024 ± 


0.017 


rv rvorv . rv ceo 

0.089± 0.553 





835± 


0.013 


6.425± 0.439 





937± 


0.015 


2.441± 


508 


NGC3938 





979± 


0.007 


1.394± 0.219 





947 ± 


0.008 


2.866± 0.250 





962± 


0.006 


1.685± 


193 


NGC4236 





730± 


0.024 


10.056± 0.779 





541± 


0.012 


16.085+ 0.402 


1 


119± 


0.027 


-3.799± 


.917 


NGC4254 


1 


052+ 


0.005 


-1.162+ 0.169 





948+ 


0.003 


2.818+ 0.110 





893+ 


0.003 


4.085+ 


117 


NGC4321 


1 


129+ 


0.004 


o 'vrv /I 1 r\ 'i A'^ 

-3.704+ 0.143 





923+ 


0.004 


3.663+ 0.142 





862+ 


0.003 


5.149+ 


113 


NGC4536 


1 


055+ 


0.009 


-1.129+ 0.289 





923+ 


0.007 


3.553+ 0.224 





924+ 


0.005 


2.963+ 


176 


NGC4559 





797+ 


0.007 


7.604+ 0.227 





910+ 


0.007 


4.109+ 0.214 





982± 


0.008 


0.971+ 


272 


NGC4569 


1 


010± 


0.010 


0.373± 0.328 





951± 


0.011 


2.692± 0.362 





911± 


0.007 


3.501± 


244 


NGC4579 


1 


103+ 


0.021 


'^/^oi . rv/^orv 

-2.681+ 0.689 





827+ 


0.011 


6.888+ 0.376 





747 ± 


0.008 


9.054+ 


264 


NGC4594 





983± 


0.009 


1.419± 0.303 





927± 


0.028 


3.644± 0.922 





888± 


0.034 


4.304± 1 


146 


NGC4625 


u 




U.Ui J 


Q 1 79_i_ n ZlQzl 







U.UiO 


3.789± 0.538 





o i Z)± 


U.Ui J 


6.698± 


500 


NGC4631 





911± 


0.004 


3.795± 0.129 





957± 


0.003 


2.568± 0.086 





980± 


0.002 


1.017±0 


075 


NGC4725 





853+ 


0.010 


5.618+ 0.348 





820+ 


0.008 


7.100+ 0.257 





849+ 


0.013 


5.577+ 


453 


NGC4736 


1 


090+ 


0.004 


-2.326+ 0.145 





998+ 


0.004 


1.188+ 0.141 





869+ 


0.002 


4.835+ 


082 


NGC4826 


1 


023+ 


0.005 


0.056+ 0.176 


1 


099+ 


0.006 


-2.197+ 0.215 





946+ 


0.003 


2.140+0 


092 


NGC5055 


1 


060+ 


0.007 


-1.395± 0.222 


1 


019+ 


0.004 


0.450± 0.141 





849+ 


0.005 


5.662± 


158 


NGC5398 


1 


214± 


0.046 


-6.073± 1.528 





697± 


0.012 


10.932± 0.386 


1 


045± 


0.017 


-1.237±0 


.569 


NGC5408 





911 + 


0.011 


4.485+ 0.366 





736± 


0.013 


9.604± 0.449 


1 


142± 


0.016 


-4.624+ 


.539 


NGC5457 





804+ 


0.006 


7.296+ 0.185 





852+ 


0.003 


5.977+ 0.113 





929+ 


0.004 


2.846+ 


136 


NGC5474 





767+ 


0.027 


8.637+ 0.896 





794± 


0.022 


7.926± 0.728 


1 


091± 


0.085 


-2.778+ 2 


.871 


NGC5713 


1 


033± 


0.007 


-0.427± 0.245 





992± 


0.008 


1.236± 0.267 





957± 


0.005 


1.789+0 


176 


NGC5866 





949+ 


0.013 


2.752+ 0.430 


1 


031+ 


0.009 


0.371+0.310 





993+ 


0.013 


0.457+ 


442 


NGC6946 


1 


011+ 


0.003 


0.285+0.104 





890+ 


0.002 


4.754+ 0.053 





903+ 


0.002 


3.728+ 


063 


NGC7331 


1 


032+ 


0.007 


-0.428+ 0.225 





985+ 


0.002 


1.625+ 0.066 





883+ 


0.004 


4.441+ 


154 


NGC7793 





851+ 


0.004 


5.742+ 0.120 





869+ 


0.004 


5.502+ 0.143 





962+ 


0.005 


1.671+0 


158 



NB: The relation is : log ^tir- a/ log S/ + b/ with Sr/^? and S/ in W kpc ^. 
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Table CI. continued 



Name 



aioo 



bioo 



ai60 



bl60 



a250 



b250 



DDO053 
Holmbergll 
IC342 
IC2574 
M81Dw 
NGC0337 
NGC0628 
NGC0855 
NGC0925 
NGC1097 
NGC1266 
NGC1291 
NGC1316 
NGC1377 
NGC1482 
NGC1512 
NGC2146 
NGC2798 
NGC2841 
NGC2915 
NGC2976 
NGC3049 
NGC3077 
NGC3184 
NGC3190 
NGC3198 
NGC3265 
NGC3351 
NGC3521 
NGC3621 
NGC3627 
NGC3773 
NGC3938 
NGC4236 
NGC4254 
NGC4321 
NGC4536 
NGC4559 
NGC4569 
NGC4579 
NGC4594 
NGC4625 
NGC4631 
NGC4725 
NGC4736 
NGC4826 
NGC5055 
NGC5398 
NGC5408 
NGC5457 
NGC5474 
NGC5713 
NGC5866 
NGC6946 
NGC7331 
NGC7793 



1.655+ 
1.615+ 
1.001 ± 
0.893± 
0.663+ 
1.038+ 
0.929± 
1.037± 
0.953± 
0.966+ 
0.946+ 
1.070+ 
1.006± 
1.053± 
0.967± 
0.963± 
0.986± 
0.985+ 
0.834± 
1.086+ 
0.966+ 
1.043± 
0.998± 
0.946± 
0.928± 
1.006± 
1.023+ 
0.975± 
0.914+ 
0.925+ 
0.985+ 
1.104± 
0.921± 
1.252+ 
0.946+ 
0.905± 
0.977+ 
0.995+ 
0.967+ 
0.947+ 
0.948± 
0.945± 
0.986± 
0.911± 
1.013± 
0.961± 
0.891 + 
1.182+ 
1.311 + 
0.994+ 
0.903+ 
0.965± 
0.985+ 
0.963± 
0.900± 
0.966+ 



0.148 
0.048 
0.003 
0.027 
0.060 
0.005 
0.004 
0.016 
0.004 
0.004 
0.015 
0.019 
0.012 
0.007 
0.006 
0.008 
0.003 
0.005 
0.006 
0.016 
0.005 
0.013 
0.005 
0.005 
0.008 
0.006 
0.009 
0.005 
0.002 
0.002 
0.002 
0.020 
0.003 
0.040 
0.002 
0.003 
0.005 
0.003 
0.005 
0.010 
0.010 
0.012 
0.002 
0.006 
0.002 
0.002 
0.002 
0.008 
0.018 
0.002 
0.019 
0.005 
0.009 
0.002 
0.002 
0.003 



-21.647+ 4.964 
-20.360+ 1.608 
0.207± 0.093 
3.851± 0.906 
11.434± 1.971 
-1.079± 0.160 
2.688± 0.128 
-1.054± 0.548 
1.862± 0.120 
1.410+ 0.131 
2.116+ 0.512 
-2.212+ 0.658 
-0.036+ 0.392 
-1.233± 0.238 
1.450± 0.195 
1.502± 0.261 
0.792± 0.098 
0.805± 0.160 
5.909+ 0.188 
-2.708+ 0.546 
1.419+ 0.177 
-1.150± 0.443 
0.324+ 0.185 
2.105± 0.161 
2.698± 0.257 
0.082± 0.211 
-0.480± 0.295 
1.097+0.177 
3.233+ 0.056 
2.816+ 0.071 
0.751+0.073 
-3.288± 0.675 
2.944± 0.114 
-8.214± 1.340 
2.134± 0.085 
3.517± 0.118 
1.099± 0.171 
0.418+0.109 
1.398+0.165 
2.023+ 0.325 
1.970± 0.335 
2.123± 0.416 
0.712± 0.053 
3.261± 0.206 
-0.235± 0.072 
1.553± 0.086 
4.024+ 0.053 
-5.840+ 0.281 
-10.214+ 0.605 
0.443+ 0.084 
3.477± 0.635 
1.494+ 0.163 
0.587+ 0.313 
1.523± 0.059 
3.711± 0.073 
1.384± 0.098 



1.749+ 
1.809+ 
1.083± 
1.235± 
1.866± 
1.058+ 
1.052± 
1.066± 
1.077± 
1.137+ 
0.952+ 
1.357+ 
1.118± 
1.074± 
l.OOli 
1.131± 
1.083+ 
1.028+ 
0.978± 
1.136+ 
1.092+ 
1.109± 
1.219+ 
1.079± 
0.979± 
1.049± 
1.061+ 
1.181± 
1.080+ 
1.078+ 
1.074+ 
1.212± 
1.029+ 
1.205± 
1.070± 
1.081 + 
1.145± 
1.005+ 
1.059+ 
1.174+ 
0.898± 
1.157± 
1.036± 
1.057± 
1.206± 
l.OlOi 
1.033+ 
1.407+ 
1.288+ 
1.044± 
0.884± 
1.066± 
0.939+ 
1.116± 
1.018± 
1.055+ 



0.190 
0.044 
0.005 
0.067 
0.266 
0.008 
0.005 
0.020 
0.004 
0.004 
0.022 
0.039 
0.026 
0.019 
0.007 
0.019 
0.005 
0.010 
0.007 
0.026 
0.008 
0.020 
0.007 
0.008 
0.012 
0.011 
0.020 
0.006 
0.002 
0.004 
0.004 
0.031 
0.004 
0.039 
0.003 
0.003 
0.008 
0.005 
0.009 
0.022 
0.017 
0.012 
0.003 
0.008 
0.003 
0.004 
0.003 
0.029 
0.025 
0.004 
0.014 
0.006 
0.008 
0.002 
0.003 
0.005 



-24.095+ 6.249 
-26.203+ 1.438 
-2.542± 0.156 
-7.464± 2.227 
-28.235± 8.779 
-1.580+ 0.274 
-1.440± 0.183 
-1.817± 0.679 
-2.278+ 0.148 
-4.363+0.142 
2.217+ 0.732 
-11.666+ 1.325 
-3.702± 0.881 
-1.516± 0.657 
0.504± 0.252 
-4.1 17± 0.648 
-2.359± 0.162 
-0.373+0.351 
0.958± 0.246 
-4.158+ 0.858 
-2.821+0.274 
-3.198± 0.681 
-7.009± 0.226 
-2.373± 0.268 
0.969± 0.419 
-1.334± 0.375 
-1.516+ 0.690 
-5.829± 0.216 
-2.435+ 0.054 
-2.351+0.122 
-2.208± 0.135 
-6.712± 1.057 
-0.657+ 0.137 
-6.446± 1.312 
-2.051± 0.104 
-2.467± 0.108 
-4.557± 0.283 
0.153+0.158 
-1.738+ 0.323 
-5.672+ 0.757 
3.687± 0.564 
-5.01 1± 0.404 
-0.889± 0.088 
-1.706± 0.268 
-6.660± 0.116 
0.020± 0.125 
-0.848+ 0.108 
-13.154+ 0.977 
-8.880+ 0.833 
-1.173+0.146 
4.251+0.477 
-1.846± 0.211 
2.340± 0.282 
-3.664± 0.074 
-0.353+ 0.114 
-1.541+0.183 



1.741+ 
1.382+ 
1.156± 
1.407+ 
0.545+ 
1.212+ 
1.212± 
1.146± 
1.238± 
1.310+ 
1.044+ 
1.657+ 
1.107± 
1.269± 
1.075± 
1.215± 
1.186± 
1.154+ 
1.130+ 
1.293+ 
1.265+ 
1.397± 
1.433+ 
1.189± 
0.994± 
1.182± 
1.212± 
1.433± 
1.266+ 
1.280+ 
1.158+ 
1.326± 
1.214+ 
1.172± 
1.218± 
1.226+ 
1.368± 
1.157+ 
1.090+ 
1.482+ 
0.849± 
1.263± 
1.102± 
1.176± 
1.436± 
1.038± 
1.206+ 
1.690+ 
1.555+ 
1.216+ 
1.006± 
1.187± 
0.979± 
1.259± 
1.160± 
1.284± 



0.118 
0.085 
0.006 
0.120 
0.214 
0.015 
0.008 
0.022 
0.008 
0.008 
0.029 
0.062 
0.038 
0.022 
0.014 
0.031 
0.012 
0.017 
0.012 
0.024 
0.011 
0.053 
0.015 
0.014 
0.025 
0.018 
0.033 
0.017 
0.003 
0.005 
0.007 
0.043 
0.008 
0.044 
0.005 
0.007 
0.014 
0.011 
0.015 
0.038 
0.020 
0.025 
0.005 
0.015 
0.007 
0.005 
0.007 
0.087 
0.033 
0.006 
0.023 
0.009 
0.008 
0.004 
0.005 
0.008 



-23.183+3.820 
-11.428+2.761 
-4.463± 0.203 
-12.589± 3.963 
15.739± 6.977 
-6.177+ 0.498 
-6.297± 0.272 
-3.888± 0.721 
-7.217± 0.280 
-9.630+ 0.267 
-0.247+ 0.985 
-20.849± 2.059 
-2.719± 1.252 
-7.336± 0.745 
-1.422± 0.456 
-6.398± 1.015 
-5.207± 0.419 
-4.021 ±0.562 
-3.745+ 0.402 
-8.743+ 0.800 
-8.138+ 0.374 
-12.287± 1.760 
-13.503± 0.491 
-5.570± 0.454 
0.926± 0.836 
-5.383± 0.595 
-5.970± 1.110 
-13.756+ 0.585 
-8.261+0.118 
-8.662+ 0.182 
-4.479+ 0.251 
-9.831± 1.418 
-6.363+ 0.264 
-4.842± 1.436 
-6.520± 0.173 
-6.855± 0.241 
-11.546± 0.456 
-4.496+ 0.368 
-2.264+ 0.505 
-15.509± 1.264 
5.719± 0.671 
-7.953± 0.844 
-2.627± 0.183 
-5.211± 0.506 
-13.737± 0.225 
-0.346± 0.158 
-6.271+ 0.241 
-21.910+ 2.893 
-16.951+ 1.076 
-6.463+0.194 
0.601± 0.755 
-5.344± 0.300 
1.527± 0.275 
-7.913± 0.124 
-4.689± 0.176 
-8.769± 0.257 
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Table C2. Calibration coefficients derived for individual galaxies to predict the TIR brightness from a combination of the 24 to 
250 fim brightnesses. 



Name C24 C70 cioo ci6o C250 



DDO053 


1 


.998+ 


0.287 


0.386+ 


0.109 


0.512+ 


0, 


,125 


0.852+ 0.264 


-0.427+ 0.412 


Holmbergll 


1 


.356+ 


0.205 


0.775+ 


0.131 


0.162+ 


0, 


,151 


0.434+ 0.328 


0.265+ 0.786 


IC342 


1 


.338+ 


0.104 


0.820+ 


0.035 


0.390+ 


0, 


,042 


0.148+ 0.059 


1.721+ 0.122 


IC2574 


1 


.847+ 


0.321 


0.486+ 


0.094 


0.531 + 


0, 


,161 


0.282+ 0.270 


0.854+ 0.604 


M81Dw 


3 


.691 + 


5.133 


0.458+ 


0.357 


0.231 + 


0, 


,207 


0.929+ 0.235 


-0.722+ 0.936 


NGC0337 


2. 


.638+ 


0.519 


0.273+ 


0.104 


0.382+ 


0. 


,191 


1.196+ 0.261 


-0.683+ 0.432 


NGC0628 


1 


.634+ 


0.219 


0.468+ 


0.068 


0.515+ 


0. 


,074 


0.508+ 0.124 


0.950+ 0.235 


NGC0855 


4 


.530+ 


1.160 


0.448+ 


0.111 


0.262+ 


0. 


,121 


0.538+ 0.188 


0.104+ 0.493 


NGC0925 


1 


.744± 


0.272 


0.389± 


0.079 


0.497± 


0, 


,087 


0.595± 0.135 


0.666± 0.210 


NGC1097 


1 


.882+ 


0.350 


0.568+ 


0.081 


0.327+ 


0, 


,115 


0.724+ 0.160 


0.498+ 0.263 


NGC1266 


3 


.419+ 


0.335 


0.388+ 


0.083 


0.295+ 


0, 


,100 


0.077+ 0.181 


0.018+ 0.551 


NGC1291 


1 


.988+ 


0.700 


0.674+ 


0.172 


0.311 + 


0, 


,163 


0.459+ 0.347 


0.563+ 0.889 


NGC1316 


1 


.457+ 


0.344 


0.792+ 


0.133 


0.278+ 


0, 


,176 


0.247+ 0.263 


1.225+ 0.560 


NGC1377 


1 


.630+ 


0.227 


0.715+ 


0.172 


0.440+ 


0, 


,212 


1.221+ 0.324 


-1.959+ 0.660 


NGC1482 


1 


.547+ 


0.274 


0.565+ 


0.091 


0.634+ 


0. 


,110 


0.816+ 0.204 


r\ A . r\ A A t 

-0.472+ 0.441 


NGC1512 


1 


.753+ 


0.282 


0.550+ 


0.087 


0.451+ 


0. 


,097 


0.336+ 0.178 


1.180+ 0.380 


NGC2146 


1 


.550+ 


0.304 


0.645+ 


0.137 


0.480+ 


0. 


,203 


0.602+ 0.171 


0.569+ 0.280 


NGC2798 


1 


.349+ 


0.257 


0.872+ 


0.193 


0.134+ 


0. 


,320 


\).151± 0.272 


0.611+ 0.515 


NGC2841 


2, 


.428+ 


1.076 


0.461 + 


0.168 


0.394+ 


0, 


,184 


0.171+ 0.307 


1 CCA, C\ A f\f\ 

1.554+ 0.499 


NGC2915 


1 


.736+ 


1.277 


0.737+ 


0.264 


0.253+ 


0, 


,257 


f\ A A 1 r\ A ^ 

0.441+ 0.465 


0.806+ 1.040 


NGC2976 


1 


.570+ 


0.443 


0.595+ 


0.134 


0.627+ 


0, 


,146 


0.040+ 0.152 


1.443+ 0.237 


NGC3049 


2, 


.253+ 


0.363 


0.478+ 


0.225 


0.005+ 


0, 


,333 


1.406+ 0.352 


-0.216+ 0.428 


NGC3077 


1 


.823+ 


0.567 


0.700+ 


0.134 


0.496+ 


0, 


,180 


0.043+ 0.267 


1.647+ 0.508 


A.Ti^/^0 10/1 

NGC3184 


1 


.464+ 


0.418 


0.718+ 


0.132 


0.375+ 


0. 


,142 


0.327+ 0.207 


A A A ^ , f\ A f\A 

1.442+ 0.401 


NGC3190 


3 


.057+ 


1.748 


0.525+ 


0.140 


0.411+ 


0. 


,186 


0.241+ 0.185 


1.066+ 0.493 


NGC3198 


1 


.777+ 


0.232 


0.485+ 


0.116 


0.431+ 


0. 


,164 


0.443+ 0.171 


1.074+ 0.220 


NGC3265 


1 


.523+ 


0.781 


0.692+ 


0.278 


0.576+ 


0. 


,407 


0.450+ 0.374 


0.018+ 0.646 


NGC3351 


1 


.984+ 


0.480 


0.553± 


0.158 


0.300± 


0. 


,176 


0.536+ 0.173 


0.976± 0.326 


NGC3521 


2, 


.042± 


0.475 


0.512± 


0.097 


0.407 ± 


0, 


,105 


0.684± 0.182 


0.529± 0.287 


NGC3621 


1 


.660+ 


0.326 


0.458+ 


0.089 


0.607+ 


0, 


,086 


0.412+ 0.148 


f\ r\f\ A . f\ ^ /" A 

0.904+ 0.264 


NGC3627 


1 


.615+ 


0.334 


0.476+ 


0.105 


0.494+ 


0, 


,118 


0.858+ 0.182 


0.004+ 0.379 


NGC3773 


2, 


.055+ 


0.673 


0.532+ 


0.224 


0.542+ 


0, 


,206 


0.094+ 0.241 


1.116+ 0.583 


NGC3938 


1 


.350+ 


0.345 


0.518+ 


0.132 


0.585+ 


0, 


,118 


0.472+ 0.207 


0.826+ 0.317 


NGC4236 


1 


.362+ 


0.136 


0.535+ 


0.074 


0.436+ 


0, 


,118 


0.633+ 0.162 


0.292+ 0.220 


NGC4254 


2. 


.389+ 


0.444 


0.551 + 


0.124 


0.515+ 


0. 


,159 


0.176+ 0.220 


1.542+ 0.368 


NGC4321 


1 


.730+ 


0.458 


0.614+ 


0.114 


0.359+ 


0. 


,126 


0.477+ 0.163 


1.167+ 0.307 


NGC4536 


1 


.665+ 


0.354 


0.488+ 


0.113 


0.608+ 


0. 


,137 


0.531+ 0.196 


0.482± 0.322 


NGC4559 


1 


.873± 


0.316 


0.396± 


0.074 


0.446± 


0, 


,101 


0.701+ 0.117 


0.502+ 0.187 


NGC4569 


1 


.093+ 


0.510 


0.935+ 


0.184 


0.389+ 


0, 


,171 


0.192+ 0.196 


1.471+ 0.414 


NGC4579 


1 


.952+ 


0.627 


0.709+ 


0.197 


0.218+ 


0, 


,219 


0.396+ 0.255 


1.256+ 0.377 


NGC4594 


1 


.582+ 


0.401 


0.750+ 


0.134 


0.172+ 


0, 


,187 


0.434+ 0.163 


1.055+ 0.159 


NGC4625 


3 


.317+ 


1.414 


0.555+ 


0.466 


0.311 + 


0, 


,541 


0.413+ 0.715 


0.930+ 1.580 


NGC4631 


2, 


.640+ 


0.358 


0.481 + 


0.060 


0.394+ 


0, 


,070 


0.719+ 0.117 


0.142+ 0.200 


NGC4725 


2, 


.688+ 


0.452 


0.633+ 


0.123 


0.128+ 


0. 


,109 


0.473+0.170 


1.195+0.275 


NGC4736 


1 


.543+ 


0.508 


0.817+ 


0.079 


0.386+ 


0. 


,102 


0.359+ 0.169 


1.064+ 0.362 


NGC4826 


2, 


.487+ 


0.465 


0.440+ 


0.139 


0.643+ 


0. 


,131 


0.085+0.130 


1.064+ 0.454 


NGC5055 


1 


.964± 


0.359 


0.477± 


0.085 


0.403± 


0. 


,103 


0.530± 0.157 


0.950± 0.231 


NGC5398 


1 


.446+ 


0.351 


0.616+ 


0.250 


0.515+ 


0, 


,348 


0.419+ 0.648 


0.644+ 1.041 


NGC5408 


1 


.699+ 


0.231 


0.808+ 


0.198 


0.270+ 


0, 


,365 


0.281+ 0.577 


0.303+ 1.607 


NGC5457 


1 


.932+ 


0.085 


0.404+ 


0.027 


0.449+ 


0, 


,043 


0.594+ 0.054 


0.755+ 0.087 


NGC5474 


1 


.793+ 


0.861 


0.503+ 


0.125 


0.312+ 


0, 


,212 


0.473+ 0.323 


1.305+ 0.668 


NGC5713 


2, 


.005+ 


0.428 


0.585+ 


0.154 


0.388+ 


0, 


,158 


0.505+ 0.302 


0.785+ 0.615 


NGC5866 


4 


.732+ 


1.320 


0.534+ 


0.159 


0.212+ 


0. 


,144 


-0.086+ 0.173 


2.137+ 0.525 


NGC6946 


1 


.478+ 


0.245 


0.695+ 


0.053 


0.456+ 


0. 


,070 


0.351+0.126 


1.124+ 0.238 


NGC7331 


3 


.376+ 


0.823 


0.389+ 


0.079 


0.396+ 


0. 


,089 


0.217+ 0.118 


1.465+0.219 


NGC7793 


1 


.601+ 


0.375 


0.537+ 


0.074 


0.467+ 


0. 


,093 


0.676+ 0.132 


0.351+0.203 



NB: The relation is : Stir= ^ c/ S/ with Stir and S/ in W kpc 
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Figure Dl. TIR brightnesses maps of the Kingfish sample in L© kpc ^ (log scale, resolution of SPIRE 250 jim) 
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Figure Dl. continued 
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Figure Dl. continued 
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